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ABSTRACT 

Ruprccht 147 is a hitherto unappreciated open cluster that holds great promise as a standard in 
fundamental stellar astrophysics. We have conducted a radial velocity survey of astrometric candidates 
with Lick, Palomar, and MMT observatories and have identified over 100 members, including 5 blue 
stragglers, 11 red giants, and 5 double-lined spectroscopic binaries (SB2s). We estimate the cluster 
metallicity from spectroscopic analysis, using Spectroscopy Made Easy (SME), and find it to be 
[M/H] = +0.07 ± 0.03. We have obtained deep CFHT/MegaCam g'r'i'z' photometry and fit Padova 
is ochrones to th e (g' — i') and 2MASS (J — Kg) CMDs, using the r 2 maximum-likelihood procedure 
of iNavlorl ([20091 ). and an alternative method using 2D cross-correlations developed in this work. We 
find best fits for isochrones at age t — 2.5 ± 0.25 Gyr, to — M = 7.35 ±0.1, and Ay = 0.25 ± 0.05, with 
additional uncertainty from the unresolved binary population and possibility of differential extinction 
across this large cluster. The inferred age is heavily dependent by our choice of stellar evolution model: 
fitting Dartmouth and PARSEC models yield age parameters of 3 Gyr and 3.25 Gyr respectively. At 
^300 pc and ~3 Gyr, Ruprecht 147 is by far the oldest nearby star cluster. 
Subject headings: open clusters: general — open clusters: individual (Ruprecht 147) 



1. INTRODUCTION 

The observational foundations of stellar astrophysics 
are studies of the Sun and stellar clusters. A few "bench- 
mark" clusters form the basis of our understanding of 
stellar evolution, and the effects of abundance, age, and 
mass on stars. When fully characterized with precise 
ages, distances and metallicities, these clusters become 
touchstones for similar stars in the field and test models 
of stellar evolution and structure. 

Galactic gravitational tidal forces are effective at dis- 
rupting most Galactic clusters on a time scale of a few 
hundred Myr (|Soderblomll2010D , so most clusters tend to 
be relatively young. This is fortunate for studies of early 
stellar evolution and massive stars: for such work stellar 
astrophysicists have access to several nearby young clus- 
ters (e.g. Pleiades ^100 - 200 Myr; Hyades and Praesepe 
~700 Myr). 

Studies of the older cool stars (age > 1 Gyr) that 
typify the field must rely on rarer and thus more dis- 
tant clusters. Studies of the typical rotation, activity 
level, and photometry of G, K and M dwarfs as a func- 
tion of age and mass, such as the WIYN Open Clus- 
ter S urvey (WOCS 7 ), the Palomar Transient Factory 
(PTF.IAgiieros et alll2011l) and the Kepler Cluster Study 
(jMeibom et alll20TTl ). investigate clusters with distances 
of 1 - 4 kpc. These larger distance moduli can make 
spectroscopic study of their low mass members extremely 
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difficult. 

Fortunately, Dias et al. (2001) and IKh archenk o et al.l 
( 2005) used catalog data to identify Ruprecht 147 (R147 
= NGC 6774), and estimated its age to be ^2.5 Gyr at a 
distance of 175 - 270 pc, making R147 by far the oldest 
nearby cluster 8 (e.g. WEBDA lists NGC 752 at 1.1 Gyr 
and 457 pc, Figured]). 

1.1. Pre-2000 literature 

Despite its promising scientific potential due to the 
unique combination of its age and distance, and despite 
having a similar distance and size to Praesepe, R147 
was complete ly overlook ed b y stellar astronom ers until 
the works by IDias et ail and IKharchenko et al.l . This is 
likely because its proximity makes R147 a very sparse 
cluster on the sky: there are only ~50 members with 
V < 11 and only ~10 with V < 9 spread over 5 square 
degrees. Its presence is also obscured by its location in 
the Galactic plane (—14° < b < —12°, in Sagittarius), 
and the fact that due to its age, it lacks the many bright 
A stars that made similarly nearby clusters so obvious, 
even to the astronomers of antiquity. 

In fact, a complete pre-2000 bibliography of R147 con- 
sists almost exclusively of entries in various catalogues. 
R147 was originally discovered in 1830 by John Herschel, 
who describe d it as "a very l arge straggling space full of 
loose stars" (IHerschell Il83l. and labelled it GC 4481 
(|HerscheWl8~63i r Since then it has appeared with numer- 
ous designations including NGC 6 774, OCL 65, Lund 883 
(|Dreverl 118881 lAlter et all H951 ILvnga fc Palousl U987t 

8 One WEBDA cluster, Loden 1, is plotted as an open cir- 
cle in Figure 1, and has properties that are apparently similar to 
Ruprecht 147. The membership and p roperties of Loden 1 were 
determined by Kharchenko ct al. (2005). They identify only nine 
1-<t members from proper motions and photometry, and none have 
measured radial velocities. The Loden 1 grouping has not been 
confirmed as a real o pen cluster, and the prop erties derived by the 
automated search of Kharchenko ct al. (2005) arc thus unreliable. 
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Figure 1. Data fro m the WEBDA database 

(Mermilliod & Paunzcn 2003) showing all known clusters in 
distance - age space. R147 is by far the oldest nearby cluster, 
and holds gr eat promise as a stand ard in fundamental stellar 
astrophysics: Kharchcnko ct al. (2005) values shown in blue, ours 
in red. WEBDA lis ts an age for M67 at 2 .5 Gyr, but we plot it at 
4 Gyr according to Pichardo ct al. (2012) and references therein. 
The open circle next to R147 denotes Loden 1, an unconfirmed 
grouping of stars with unreliable properties - see footnote 8 for a 
discussion. 



19*18 m 19 h 16 m 19 h ' 

Right Ascension 



Figure 2. This astrophotograph of a portion of the Ruprecht 147 
field was taken and kindly provided by Chris Beckett and Stefano 
Mcncguolo of the Royal Astronomical Society of Canada. We have 
attached an approximate coordinate system solution (we have not 
solved for the field distortions) and circled the 80 high-confidence 
members in red. There are 52 stars in this image with V < 10, and 
only 17 are members of R147 and are circled green. Of the 47 stars 
with V < 9 within ?s 2° of the cluster center (extending beyond 
this image), only 11 are members. It is remarkable that Hcrschcl 
correctly identified this as an open cluster in 1830. 



IMermilhodl I1995D . Some star charts have even des- 
ignated R147 as an asterism, and not a true cluster 
(e.g. "Burnham's Celestial Handbook: An Observer's 
Guide to the Universe Beyond the Solar Sys tem" lists 
NGC 6774 as "possibly not a true cluster" (B urnhaml 
1966)). The name we use here originates from iRuprechtJ 
(1966), who classifi ed R147 as a II I-2-m cluster in the 
Tr umpler system dTrumplerl 119301 160). According 
to lArchinal fc Hvnesl (|2003|) T Brian Skiff realized that 
NGC 6774 and R 1 47 ar e likely the same star cluster. 
lArchinal fc Hvnesl (|2003l ) describe R147 as a "45' sized 
V-shaped group of bright stars" that is "a sparse pos- 
sible open cluster", and estimate the cluster center as 
the location of HD 180228 (while this star's photometry 
apparently places it on the R147 red giant branch, the 
Tycho-2 proper motions, -1.6 and -6.3 mas/yr in right 
ascension and declination, are inconsistent with cluster 
membership, see Figure [3]). Figure [2] highlights our high- 
confidence members on an optical image. Herschel's clus- 
ter identification is truly amazing, given the lack of a well 
defined cluster core. But those arguing for the asterism 
status were not entirely wrong either: of the 51 stars 
with V < 9 within m 2° of the cluster center, we confirm 
only 11 as members. 

1.2. Recent work in the literature 

Only in the last decade has R147 receive d any indi- 
vidual attention in studies of open clusters. [Dias et all 
(|2001l) first identified R147's membership based on the 
stellar population's common pro per motion : selecting 
stars in the Tycho-2 Catalogue (jHog et all 12000) that 
were spatially coin cident with BDA clusters (The Open 
Cluster DataBase, Mermilliod 1995), they determined 
cluster membership with the Tycho-2 proper motions us- 
ing the statistical method of iSandersi (|197lD and found 
33 stars with mean prope r motion of M a = — 0.8±2.3 and 
(is = — 28.5±2.3 mas/yr. iDias et ail (|2001l) also provided 
the first distance estimate based on only two Hippar- 



cos p arallax measurements 9 (HIP1, iPerrvman fc ESAl 
[19971) : 7T = 3.57 ± 1.01 mas (280 ± 79 pc) for HIP 94635 
(CWW 1) 10 , and 3.75 ± 1.04 mas (267 ± 74 pc) for HIP 
94803 (CWW 2), which they average to 3.7 ± 0.2 mas, 
estim ating the distance to R147 to be 250 pc 11 . Since 
then, Ivan Leeuwenl (|2007af) (HIP2) has performed a new 
data reduction and issued an updated catalog with par- 
allaxes of 5.48 ± 0.65 mas (182 ± 22 pc) for HIP 94635, 
a nd 4.92 ± . 79 ma s (203 ± 33 pc) for HIP 94803 12 . 

IDias et all (|2002l ) compiled all available data for 2095 
galactic clusters ( The New Catalogue of Optically Visible 
Open Clusters and Candidates, or DAML02) and pub- 
lished an updated membership list and cluster properties 
for R147: 25 members, proper motion fi a = —0.9 ± 0.3 
and /is = —29.3 ± 0.3 mas/yr, RV = 41 km s" 1 ( from 
the single published measurement in iWilsonl ()1953l ). see 
t }3~2|) . distance = 200 pc, color excess E(B - V) = 0.2 
mag., and an age 3.2 Myr (presumably from misiden- 
tifying blu e stragglers as main sequence turnoff stars). 
IDias et ail re-classified R147 as IV-2-p (Trumpler sys- 
tem) . 

Following their 2002 work, IDias et al.l (|2006l) se- 
lected all clusters in their DAML02 catalogue with 
known distances and queried the UCAC-2 catalogue 
(jZacharias et al.l l2004bf ) for all stars within the mea- 

9 Actu ally, three R147 members appear in the Hipparcos catalog, 
see 34~7l 

10 Throughout this paper, we will refer to individual stars 
with the designation "CWW #" (CWW = Curtis, Wolfgang and 
Wright). Our membership list provides 2MASS IDs, astrometry, 
photometry, radial velocities, and membership probabilities for 108 
stars. The CWW ID numbers sort these stars according to V mag- 
nitude (see Table 6). 

11 although this is a numerical error as 1000 / 3.7 = 270, not 
250. 

12 Ivan Lccuwcn 1 20 07bl . §3.3.1, 3.3.2) cautions against deriving 
distances and distance moduli from parallaxes when the relative 
error is greater than 10%. The Lutz-Kelker bias can also introduce 
a 0.1 magnitude systematic offset at 10% relative error 
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sured cluster radii, plus 2', of their tabu lated clus- 
ter ce nters. Employing similar methods as Di as et ahl 
(2001), they derived a mean proper motion for R147 of 
[i a — —4.6 ± 0.4 and fis — —5.6 ± 0.4, and identified 
200 cluster members. Figure [3] shows the proper motions 
for stars in the R147 field , color-shaded by membership 
probability as derived by iDias et al.l (|2006l ). The black 
cir cle highlights the prope r motion of R147 according 
to Kharchenko et al. (2005) and confirmed in this work, 
and shows that the Dias algorithm missed the c luster , 
locating the field stars instead. The IDias et all ([20061 ) 
mem bership list and c luster parameters are thus unreli- 
able. IDias et ahl (|2006l ) attribute their algorithm's failure 
to the large angular size of R147. 
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Figure 3. Proper motion diagram of stars in the R147 field, 
color shaded by membership probability as derived by IDias et alJ 
(200 6). The black circ le highlights the proper motion of R147 
I Kharc henko ct al. 2005], and confirmed here). The Dias member- 
ship probablities are clearly in error. 



A similar autom ated effort has been undertaken by 
iKharchenkol (|2001| ). who assembled the All-Sky Com- 
piled Catalogue of 2.5 Million Stars (ASCC-2 5), includ- 
ing p roper motions from the Tycho-2 catalog (lH0g et ahl 
2000), Johnson BV photometry, and radial velocities and 
s pectral types when t hey a re available. 

iKharchenko et all (|2005l ) searched this catalog and 
identified 520 Galactic open clusters, including R147. 
Their algorithm determined the core and cluster an- 
gular radii, and the distances, mean space motions 
(proper motion and radial velocity), and ages of the 
clusters. Three i mportant differences exist between the 
IDias et all (I2006D members hip a nd properties and t hose 
of IKharchenko et a l. (2005): (1) IKharchenko et al.l cor- 
rectly identify the cluster, cataloging 41 1-a members; 
(2) they provide the first reliable age estimate of 2.45 
Gyr from their isochrone fitting; and (3) they claim a 
new distance of only 175 pc, 75 pc closer than that in- 
ferred from the original Hipparcos parallaxes, but sim- 
ilar to the distances derived in HIP2. While we deter- 
mine a similar age of ~2.5 Gyr, we derive a distance 
d w 300 pc (i )4.31 I4.4[) by fitting isochrones to a spec- 
troscopically derived T e g - logg diagram, and 2MASS 
(J - K s ) and CFHT/MegaCam (g' - i') color - mag- 
nitude diagrams. Figu re [4] plots the CMD used by 
IKharchenko et all (|2005|) to derive age and distance. The 



Tycho-2 BV photometry is magnitude limited at V ~ 11, 
near the R147 main sequence turnoff. ASCC-2.5 is sup- 
plemented with various ground based photometry for 
fainter magnitudes, which Figure|3]demonstrates is insuf- 
ficient for main sequence fitting. While the MSTO pro- 
vides a strong constraint on the age, t he descrepancy be- 
tween our derived distance and that of lKharchenko et al.l 
can be explained by the ill-defined (B — V) main se- 
quence. Their analysis was also hindered by a lack of a 
spectroscopically determined composition, and they as- 
sumed Solar metallicity. Without an accurate metallic- 
ity, and with a main sequence dominated by photometric 
error, it is difficult to disentangle vis ual extinction, age 
compo sition and distance. Instead, IKharchenko et al 
2005) assumed A v = 0.465 from the iSchlegel et al 
1998) dust map at their location for the cluster cen- 
ter, even though according to this dust map Ay varies 
from 0. 3 to 0.6 mag, ac ross the cluster (see iJ4.ll) . Al- 
though |Dia^^3jll2^1 were t he first to determine the 
distance, Schilbach etaLl (2006) were the first to discuss 
Ruprecht 147 specifically as an old nearby cluster in a 
peer-reviewed publication. 

Despite these i ssues, the works of IDias et all and 
IKharc henko et al,s are significant because they essen- 
tially re-discoverd Ruprecht 147 and provided the first 
good evidence that R147 is in fact the oldest nearby star 
cluster. 

Most recently, iPakhomov et al.l (|2009D observed three 
cluster red giants and spectroscopically measured radial 
velocities and stellar properties (discussed in iJ3.2[ 14. 2p . 
They determined the cluster metallicity to be super- 
Solar, thereby decreasing the estim ated age, from a fit 
to an enriched Pa dova isochrone (|Girardi et al.l 120001 : 
iMarigo et all 120081 ) 13 to -1.25 Gyr, and derived a dis- 
tance of 280 ± 100 pc, along with a color excess of 
E(B -V) = 0.11 (or A v = 0.34, assuming R v = 3.1). 

Ruprecht 147 has also a ppeared in the open clus ter lu- 
minosity function study of lElsanhourv et all (|2011D and a 
paper on Galacti c kinematic s and structure as defined by 
open clusters by iZhul (|2009f ). but these works undoubt- 
edly suffer from a poorly determined membership and 
uncertain cluster properties. 

We have begun an observational campagin to charac- 
terize R147, catalog its members, and prove its bench- 
mark status. Here we present our initial efforts, detail- 
ing in particular our R147 membership search that more 
than doubles the number of known cluster members (*J3J), 
and our derivation of the cluster's age, distance, and 
metallicity (Q. We begin with an overview of our pho- 
tometric and spectroscopic datasets. 

2. OBSERVATIONAL DATASETS 

Cluster members are identified by their common space 
motion, determined from proper motions and radial ve- 
locities, and by their placement on a color - magnitude 
diagram (CMD). We utilize the NOMAD, UCAC-3 and 
PPMXL astrometric catalogs for proper motions. We 
have high resolution, single order echelle spectra from 

13 The Padova isochrones are available at 
|http: //stev. oapd.inaf.it/cgi-bin/cmd We primarily use these 
stellar evolution models because the Padova group provides 
synthetic photometry in a large number of systems, including the 
CFHT/MegaCam g'r'i'z' filter set allowing us to analyze our 
optical photometry. 
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Figure 4. Ruprecht 147 color — magnitude diagrams. The left panel shows the (B — V) photometry used by Kharchcnko et al. (2005) 
to estimate age and distance by isochrone fitting. These data are magnitude limited at the main sequence turnoff. The central and right 
panels plot the same stars in black, along with our additional members in red. The main s equence is better d enned in the 2MASS NIR 
(J — Kg) and our optical (g' — i') CMDs, which explains the «80% discrepancy between the [Kharchcnko ct al. distance of 175 pc and our 
value of Rs300 pc. The g' and i' error bars are set at 0.03 magnitude. The color errors are the magnitude errors added in quadrature. 
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MMT Observatory; and high-resolution, cross-dispersed 
echelle spectra from Lick, Palomar, and Keck Obser- 
vatories. We acquired deep g'r'i' ' z' photometry of a 4 
square degree field with CFHT/MegaCam, and utilize 
NIR JHK S photometry from the 2MASS Point Source 
Catalog. Other observing projects are underway, includ- 
ing deep NIR imaging with UKIRT / WFC AM (PI Adam 
Kraus), a 250 ks exposure of the cluster core with Chan- 
dra/ACIS (PI Steve Saar), and an RV survey for K and M 
dwarf members with Magellan/MIKE+MagE (PI Steve 
Saar). 

2.1. Astrometric catalogs 

Our initial list of candidate members was drawn 
from the NOMAD and UCAC-3 catalogs. The Naval 
Observatory Merged Astrometric Dataset (NOMAD, 
iZacharias et al.lfeOtMal ) combines data (positions, proper 
motions, and BVR/ JHK ph otometry) for over 1 bi l- 
lion stars from the Hippar cos ()Perrvman fc ESAl 119971) . 
Tycho-2 (|Hpg et all [2000h. UCAC-2 (IZacharias et al.l 
I2004bl). USNO-B1.0 Monet et al.l 120031). and 2MASS 
(jSkrutskie et^l200fih catalog s. The Third USNO C CD 
Astrograph Catalog (UCAC-3. IZacharias et al~ll2009t ) ex- 
pands on NOMAD by improving UCAC-2 in many ways, 
including complete sky coverage, reduced systematic er- 
rors for CCD observations, deeper photometry (R w 8 
- 16) for ~80 million stars, and improved astrometry 
(resolved double stars, inclusion of several new catalogs, 
and re-reduction of early epoch photographic plates to 
derive proper motions). 

In this paper, we use prop er motions from the PPMXL 
catalog (jRoeser et al.ll2010l) . and provide these values in 
Table 6. PPMXL utilizes astrometry from the USNO- 
B1.0 and 2MASS catalogs to calculate proper motions in 
the ICRS system for approximately 900 million objects, 
including ~410 million with 2MASS photometry. The 
catalog covers the entire sky down to V ~ 20. PPMXL 
was released in 2010, after we had derived our initial 
membership catalog. Some stars have NOMAD and/or 
UCAC-3 proper motions consistent with cluster member- 
ship, but are discrepant according to the PPMXL values 
(and vice versa). We include these stars in our member- 
ship list despite this, and we evaluate their probability 
of membership based on all available kinematic and pho- 
tometric data (f}3]). 

2.2. Lick 3-m and Palomar 200-in spectra and radial 
velocities 

We performed initial radial velocity confirmation of 
suspected members to verify the existence of the clus- 
ter with the Hamilton echelle spectrometer on the 120 
inch Shane telescope at Lick Observatory [R ~ 50, 000; 
lVogti ll987). Our objectives were to obtain RVs of known 
and suspected members, to identify new members, and to 
obtain high resolution spectra of the brightest members 
at high signal-to-noise ratios (SNRs) for more detailed 
analysis of abundances and chromospheric activity. 

We observed candidate cluster members on UT 2007 
July 31 - August 1 and 200 7 August 22 - 23 includ- 
ing the members identified by iKharchcn ko et al.l {2005). 
To locate additional candidate members, we selected 
stars from the NOMAD catalog that were within 1°.25 
of the published cluster center and that had UCAC- 
2 and TYCHO-2 proper motions within 5 mas/yr of 



the IKharchcn ko et all value. Although there are over 
750,000 NOMAD stars in the field due to its large size 
and location in the Galactic plane, the cluster is sepa- 
rated enough from the field in proper motion space that 
this yielded a list of 1348 stars, illustrated in Figures [5] 




Figure 5. Proper motion cuts made to NOMAD stars to identify 
candidate R147 members. Of the over 750,000 NOMAD stars in the 
1°.5 radius region centered on R147, 38,623 have < J — Ks < 1 
and 4 < J < 14 which are shown in gray. There are 1348 stars with 
NOMAD proper motions within 5 mas/yr of the cluster value, 280 
of which have J < 14 and are plotted in blue. Our final membership 
list of 108 stars are circled in red. This plot demonstrates that there 
is really only one obvious sequence at bright magnitudes where 
the proper motions are reliable, including a well-defined red giant 
branch. 



To further vet this list, we used NOMAD (B ~ V) and 
2MASS (J — Kg) color - magnitude diagrams to iden- 
tify stars consistent with an assumed distance of 230 pc, 
a compromise between the HIP 1 distance to the cluster 
(270 pc) and the value of iKharchenko et al.l (175 pc). 
We combined the CMDs and proper motion informa- 
tion to estimate crude membership probabilities based 
on the Hipparcos main sequence with no reddening cor- 
rections, calculated generously to account for uncertain- 
ties in the cluster parameters and for the poor quality of 
some of the NOMAD proper motion entries, and we fa- 
vored brighter targets to improve the efficiency of vetting 
candidate members at the telescope. 

We drew from this list, sorted by membership proba- 
bility, to choose targets for spectroscopic study at Lick 
Observatory. We used these spectra to measure radial 
velocities for the stars and determine the space motion 
of the cluster. 

2.2.1. Data acquisition and raw reduction 

We adopted the spectrograph setup procedure of the 
California and Carnegie Planet Search, which placed 
bright emission lines from a thorium-argon (ThAr) lamp 
on specific pixels to approximately reproduce a known 
wavelength solution to a fraction of a pixel. Our prior 
experience using the Hamilton spectrograph at coude fo- 
cus revealed that the wavelength solution is reliable to 
a pixel or two over the course of the night. This was 
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sufficient for our purposes of measuring radial velocities 
to < 5 km s -1 , a precision which allows most interloping 
field stars to be identified, so we did not attempt any 
further wavelength calibration throughout the night. In 
practice, our radial velocity accuracy proved to be much 
better than 5 km s~ . 

Observing conditions were good, and we obtained sev- 
eral high SNR spectra of radial velocity standard stars 
of various spectral types through o ut th e nights, chosen 
from the catalog of iNidever et al.l (|2002T ) . For candidate 
members we used exposure times of 60 - 90 s, depending 
on the magnitude of the star. For fainter stars we ob- 
tained SNRs as low as 1 per pixel, which is sufficient for 
our velocity work because of the broad spectral coverage 
of the Hamilton spectrograph. This strategy allowed a 
large number of stars to be observed in our allocated time 
for this low-declination cluster, which only spent a few 
hours per night at sufficiently low airmass to be useful. 

The raw spectra were processed with the standard 
Hamilton Spectrograph data reduction pipeline used 
for precise radial velocity work by the California and 
Carnegie Planet Search, which includes bias subtraction 
and flat fielding of each frame and which results in a 
one-dimensional spectrum for each of 92 orders. 

We calculated an empirical blaze function for each or- 
der by fitting a polynomial to the spectra of several 
rapidly rotating B stars that we observed for this pur- 
pose. These stars show no high resolution spectral fea- 
tures, and we corrected orders contaminated by the ef- 
fects of the broad Balmer lines by averaging the polyno- 
mial blaze function of the neighboring orders. Variations 
in slit illumination from target to target created appar- 
ent continuum variations that were not perfectly removed 
by this process, and the nature of the polynomial fitting 
process caused the fit to diverge from the actual spec- 
trum significantly at the edges of orders. The resulting 
spectra were nonetheless sufficiently flat that the cross 
correlations required for our data analysis ( fl2.2.3j) could 
be confidently performed. 

2.2.2. Palomar spectra 

We followed a similar procedure at Palomar Observa- 
tory to determine membership probabilities and activ- 
ity measurements of fainter candidate members. Tar- 
get stars were drawn from the same sorted list that was 
compiled for the Lick observing run the previous year, 
including 25 targets that were chosen for follow-up ob- 
servations based on qualitative examination of the radial 
velocity measurements derived from the Lick data, ei- 
ther because the signal-to-noise ratio of the Lick data 
was too low for a definitive velocity measurement or the 
star showed evidence of binarity, necessitating a second 
epoch. 

We observed on 2008 Au gust 5 and 18 with the East- 
Arm Echelle (R ~ 33. 000. iPerl I1995D on the Hale 200- 
inch at Palomar, following our earlier procedure of short 
integrations at very low signal-to-noise ratio (the addi- 
tional aperture of the 200- in over the 3-m, somewhat mit- 
igated by the low throughput of the East Arm Echelle, 
allowed us to explore fainter targets, or brighter targets 
at better SNR). 

2.2.3. Radial velocity determination 



Although we adopted the rough wavelength calibration 
used for planet search work, we did not attempt to use 
this calibration to measure our radial velocities. Rather, 
we extracted radial velocities in pixel space by cross- 
correlating the spectra of our candidate cluster members 
with those of our observed RV standard stars (a sim- 
ila r pixel-space cross- correlation method was employed 
bv lNorris etaTl 1)201 If ) 1. To reduce the errors introduced 
by comparing two stars of different spectral types, we 
paired each candidate member to an RV standard star 
that minimized the difference between their V — J colors 
(A(V - J)), with V - J = 0.8 for the bluest standard 
star and V — J = 2.4 for the reddest. 

Imperfect flat fielding produced a sharp spike at ex- 
actly zero shift in the cross-correlation functions (CCFs), 
and the presence of telluric lines created a narrow peak 
there, complicating the radial velocity measurements de- 
rived from these CCFs. This justified our use of velocity 
standard stars as cross-correlation templates rather than 
high SNR spectra of actual cluster members, since the 
standards have different radial velocities than the cluster 
and so the true CCF peak is far removed from the spuri- 
ous peaks at zero shift. To further address the problem 
of telluric lines, we empirically rejected those portions 
of the spectrum where these lines dominated the CCF: 
after dividing each of the 92 orders into three segments, 
we discarded from all spectra those segments that showed 
strong telluric peaks near zero shift. 

Computing these CCFs for the different combinations 
of RV standard stars with A(V — J) < 0.5 allowed us 
to calibrate the conversion from pixel space shifts to ra- 
dial velocities, after applying a barycentric radial velocity 
correction. This calibration step thus obviates the need 
for a transformation into wavelength space. Specifically, 
we fit a linear function with zero intercept to the mea- 
sured RV standard stars' CCF pixel shifts as a function 
of the difference in their radial velocities, giving us the 
velocity shift per pixel in each spectrum segment. The 
root mean squared of the residuals of this fit is less than 
0.6 km s _ ; this provides our best measure of the sys- 
tematic velocity precision we expect at high SNRs. 

Comparison of this calibration constant among the seg- 
ments confirmed that the velocity shift per pixel of the 
Hamilton spectrograph is nearly constant for each of the 
92 orders. This is not surprising given that both the reso- 
lution, R (— A/AA per resolution element), and the sam- 
pling, s (pixels/resolution element), are nearly constant 
across an echellogram, and that our calibration constant, 
having units of km s _1 /pixel, is essentially c/Rs, where 
c is the speed of light. This allowed us to add the CCFs 
of the remaining segments together to improve the SNR 
of the stellar signal, and enabled the clear identification 
of a peak and its associated pixel shift in the combined 
CCF. To be conservative, we divided each spectrum into 
three sets of segments, corresponding to the left, middle 
and right sides of each order. After separately summing 
the CCFs in each set, we required that the location of 
the tallest peak in the summed CCFs to be identical in 
all three sets; however, when the side segments produced 
noisy CCFs, as was the case for stars with SNR ~ 10, we 
used the location of the tallest peak in the middle seg- 
ments' summed CCF, as long as the peak met our high 
quality classification. 
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We visually inspected each CCF produced by our pro- 
cedure and classified the stars into two categories based 
on the quality of their CCFs: either the combined CCF 
had one clear peak, which corresponded to high SNR 
spectra, or the combined CCF had multiple peaks of ap- 
proximately equal height, which indicated that the CCF 
was dominated by noise, and which usually corresponded 
to spectra with SNR < 2 per pixel (some of these dis- 
carded stars were later revisited in the Palomar observing 
run in order to acquire higher SNR spectra). 

We found a clear clustering of barycentric velocities 
near 43 km s _1 , which is w i thin 2 km s _1 of the clus- 
ter RV quoted bv lDias et al.l (|2002f) that was base d on a 
single measurement of a single putative member ( §3.2p . 
Upon closer inspection of the Lick velocities for signs of 
systematics, we found that these apparent cluster mem- 
bers' velocities exhibited a slight correlation with time 
from the beginning of the observing run to the end, with 
magnitude 2-4 km s _1 . We fit this trend to a lin- 
ear function and removed it. We have also observed 49 
of these stars with MMT/Hectochelle (with very high 
signal-to-noise ratio, we expect all to have RV precision 
«0.5km s~\ and measure differences between the 

Lick/Palomar and Hectochelle RVs as large as 5 km s~ , 
except for a possible single-lined spectroscopic binary 
(SB1) with a difference of 20 km s _1 . We interpret this 
offset between telescopes and spectrographs as a measure 
of systematic error in our absolute barycentric radial ve- 
locities. 

Figure [5] illustrates the resulting RVs and shows a clear 
clustering around the cluster velocity. We tentatively 
identified as cluster members any stars with a measured 
radial velocity between 32 and 54 km s~\ or roughly 
twice the typical systematic error. 

2.3. Spectra from Keck/HIRES 

Spectra of four cluster members were obtained on 2008 
September 12 and 18 and for two members on 2011 Oc- 
tober 17 with the High R esolution Echelle Spectrometer 
(HIRES. IVogt et al.lll99l on the 10-m telescope at Keck 
Observatory. The stars were kindly observed by the Cal- 
ifornia Planet Survey team (CPS ) without an iodine 
cell, and with the B5 decker (slit of 3.5" length and 
0.861" width), giving a typical resolution R ~ 50,000 
in the 3360 - 8100 A bandpass. Exposure times were 
monitored with a photomultiplier tube exposure meter to 
ensure high signal-to- noise (S/N ~ 50 — 100). These ob- 
servations are summarized in Tabled] and were re duced 
by the standard CPS pipeline. iChubak et al.l (|2012fl mea- 
sured absolute radial velocities (results dis cusse d in i )3.2l) 
and we derive stellar properties in Section 14.21 

2.4. Spectra from MMT/Hectochelle 

We obtained high-resolution spectra with 
MMT/Hectochelle in the vicinity of the Ca n H & K 
lines for 48 members (as determined from Lick/Palomar 
RVs), 10 candidate members (from astrometry and 
photometry alone), and 23 potential astrometric ref- 
erence stars. These data provide radi al velocities, 
chromospheric activity indicators (e.g. (Wri ght et al.l 

14 http:/ /exoplanets.org/cps.html 



120041) ) and gravity diagnostics (via the Wilson-Bappu 
effect, iWilson fc Vainu Bappul I1957D useful for iden- 
tifying background giants as astrometric references. 
The remaining fibers not allocated for sky subtraction 
were assigned to proper motion candidates with incon- 
sistent photometry to assess how R147 members are 
kinematically distinct from the field. 

MMT is a 6.5-m telescope located at the Fred 
Lawrence Whipple O bservatory on Mt. Hopkins, AZ 
(jFabricant et al.l 120041 ) . Hectochelle is a high- resolution 
(R ~ 32,000 - 40,000) fiber-feb spectrograph, which 
provides simultaneous observ ations for 240 targets in 
a one square de gree field (|Szentgvorgyi et al.l 11998; 
iFiiresz et al.l 120 08). We observed the central square de- 
gree with the 'Ca41' Ca II H & K filter with lxl on-chip 
binning. Eight total hours were obtained to ensure suffi- 
cient signal-to-noise for a future chromospheric activity 
study. All observed targets have g' — 9 — 15.5. 

Twelve 40-minute exposures were obtained over the 
nights of UT 2010 July 4-5. These data were reduced 
with an IRAF 15 -based automated pipeline developed at 
the Harvard-Smithsonian Center for Astrophysics, pro- 
vided and run by Gabor Fiiresz and Andrew Szentgy- 
orgyi, which flat-fielded, cosmic-ray removed, and wave- 
length calibrated our targets and sky flats. The wave- 
length solution was determined from Thorium- Argon 
(ThAr) lamp comparison spectra, with an RMS precison 
of 0.2 — 0.5 km s _1 (for reduction details, see lMink et al.1 
120071) . 

Radial velocities were measured by cross-correlating 
the target spectrum with Solar spectra obtained from 
the sky flat exposures, then corrected for Earth's helio- 
centric motion. We checked the fiber-to-fiber and day- 
to-day stability of the spectrograph by measuring veloc- 
ity shifts determined by cross-correlating matched ThAr, 
Solar, and target spectra. The fiber-to-fiber velocity shift 
on the first night was 12 ± 35 m s _1 (for the ThAr spec- 
tra) and 200 ±200 m s _1 (for the Solar spectra). We also 
measure a night-to-night variation between each fiber of 
31 ± 41 m s" 1 (ThAr) and 200 ± 180 m s'^Solar). The 
RVs measured each day for R147 stars show a mean abso- 
lute difference of 0.23 km s _ . In summary, fiber-to-fiber 
and day-to-day offsets and errors are well under or com- 
parable to the precision set by the wavelength solution. 
The RV distribution for 49 member stars is shown in 
Figure El and is discussed in £13.21 

2.5. Preliminary optical photometry 

We imaged a 4 square degree field in the optical g'r'i'z 1 
bands in 2008 A pril and May with CFHT/MegaCam 
(jHora et al J 11991 16 . With MegaCam's one square de- 
gree field of view, four fields were required to cover the 
majority of the known cluster. The fields are outlined 
over a 2MASS J band mosaic image in Figure [B] Six ad- 
ditional surrounding fields were imaged solely in i' band, 
for the purpose of first-epoch astrometry for the entire 
cluster, including any extended halo. We obtained both 
a series of short (1 second) and long (few minutes) ex- 

15 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 

1B http:/ /www. cfht.hawaii.edu/Instrumcnts/Imaging/McgaPrime/ 



Table 1 

Keck/HIRES observations. 



CWW ID 2MASS ID Obs Date Exposure Time Airmass V a S/N b Notes c 







JD 


seconds 




mag. 






72 


19165800-1614277 


2008-09-12 


210 


1.43 


11.52 


50 


G dwarf / SB2 


78 


19160879-1524279 


2008-09-12 


170 


1.44 


11.82 


60 


late F dwarf 


21 


19132220-1645096 


2008-09-18 


90 


2.17 


9.98 


80 


Subgiant 


22 


19172382-1612488 


2008-09-18 


93 


2.12 


10.04 


70 


mid-F MSTO / SB1 


41 


19164495-1717074 


2011-10-17 


167 


1.39 


10.61 


80 


mid-F dwarf MSTO 


91 


19164725-1604093 


2011-10-17 


822 


1.32 


12.39 


50 


early G dwarf 



-15°00' 



-15°30' 



-16°00' 



a V magnitudes arc drawn from the NOMAD catalog. 

Signal-to-noise ratio measured in the spectral order encompassing the Mg b triplet, in the 5034 — 5036 A continuum. 
c Approximate classification, performed by matching spectroscopic and photometric properties to isochronc masses ( S|3.3t ■ 
SB1 status suggested by inconsistent RVs from multiple epochs; MSTO — main sequence turnoff 

served in all the filters utilized on a given night, but the 
zero point scatter across an observing run for each filter 
ranged from 0.0073 to 0.0180, and is therefore quite sta- 
ble. Frames obtained on photometric nights provide the 
means to calibrate observations taken under less trans- 
parent conditions (the image scaling is done by TER- 
APIX, see below) . The zero points were determeind after 
the application of the superflat, and therefore are valid 
for all 36 CCDs. The photometric and astrometric cali- 
bration data are stored in the FITS image headers, and 
the data transferred to the Canadian Astronomy Data 
Centre (CADC) in Victoria 18 . 

TERAPIX performed the final photometric and astro- 
metric reductio n of our MegaCam imagery 19 . The TER- 
APIX pipeline (|Bertin et al.l l2002t) takes the detrended 
images and the preliminary calibration from Elixir, and 
completes a final photometric and astrometric calibration 
and provides source merged catalogs. First the images 
are re-scaled: the photometry is analyzed in each over- 
lapping frame, and the frames are re-scaled to the pho- 
tometry in the image with highest flux per source, which 
is considered to be the least extinguished. The overlap- 
ping images are stacked and the sources are re-extracted. 
The TERAPIX pipeline co-addition and astrometric cal- 
ibration modules are now maintained at AstrOmatic.net, 
and document ation for each can be found at the SWarp 20 
and SCAMP (|Bertinl 120061 Software for Calibrating As- 
troMetry and Photometry,) 21 webpages. 

TERAPIX handles the CFHT Legacy Survey reduc- 
tion. The CFHTLS 22 and TERAPIX's CFHTLS reduc- 
tion 23 webpages provide additional details relevant to the 
final photometric and astrometric calibrations. 




ig^tr i9"ia m ig^e" 1 i9'u" 

Right Ascension 



19"12 m 19 h 10 m 



Figure 6. The four fields we imaged in g'r'i'z' with 
CFHT/MegaCam, overlaid on a 2MASS J band mosaic image gen- 
erated with Montage (http://montage.ipac.caltech.edu). We log- 
scaled and smoothed the image with a 3" boxcar. 



posures in queued service observing mode. Typically, 5 
dithered exposures were obtained for each field and ex- 
posure time. 

These observati ons were pre-processed at CF HT with 
the Elixir pipeline (jMagnier fc Cuillandrell2002fl 17 . Elixir 
creates master bias, dark and flat images, which are used 
to detrend the observation frames. SExtractor identifies 
sources and determines their pixel coordinates and raw 
flux. The astrometric calibration is performed by com- 
parison to the USNO-B1.0 catalog. 

Photometric magnitudes are calibrated from the in- 
strumental magnitude with the application of a zero 
point, an airmass term and a color term. The coefficients 
are derived from observations of standard stars. Every 
night, one Landolt (1992) field was observed (SA-101, 
SA-107 and SA-113), along with two spectrophotomet- 
ric standards (i.e. an O star or white dwarf: Feige 110, 
CD 153, HZ 43, BD+28 4211) and at least one CFHTLS 
Deep field. The zero points for each frame were deter- 
mined from 13 to 26 standards observed in 4 to 9 separate 
images during a run. Standards were not necessarily ob- 

17 http:/ /www. cfht.hawaii.edu/Instruments/Elixir/home. html 



TERAPIX kindly provided us merged g'r'i'z' source 
catalogs for each field and exposure duration. The short 
exposures saturate at g' w 9.5 and the long exposures 
saturate at g' ~ 16. Sources with g' ~ 17 — 18 have 
consistent photometric magnitudes in each catalog, and 
sources are detected down to g' ~ 24 in the long exposure 
catalog. We therefore have photometry covering 10 < 
<?' < 24. _ 

Our faintest red giant branch member is g' = 10.53 

18 http:/ /wwwl.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/cadc/ 

19 TERAPIX is a data reducti on center located at the Institut 
d'Astrophysique in Paris, France: http://terapix.iap.fr/ 

20 http: / / www.astromatic.net / soft ware/s warp 

21 http: //www. astromatic.net / software/scamp 

22 http://www.cfht.hawaii.edu/Science/CFHTLS/ 

23 http:/ /terapix.iap.fr/article.php?id_article=383 
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and i' = 9.73. Given the saturation limit in each band at 
approximately 9.5, the majority of the red giant branch 
stars are saturated in g' and i'. We will include the 
optical RGB in our figures in this work, but with stars 
plotted with open circles to distinguish these data from 
the more reliable optical photometry for the rest of the 
membership. 

We estimate the photometric error by making use of 
the overlapping regions between the four imaged fields 
(see Figure [S]). We matched all stars in the overlap re- 
gions in our bright source catalog (short exposures), and 
find a total of 1575 unique sources with g' < 18. Figure[7] 
plots the mean versus standard deviation of g' for the 2 - 
4 independent measurements, depending on the number 
of overlapping regions containing the source. We find a 
typical value of A g > = 0.035 mag., but this is probably 
larger than what should be assumed for the photometric 
precision across the field, because one of the sources usu- 
ally lies very close to the edge of one of the fields, where 
the mosaic dithering is incomplete and the photometry 
less reliable. 

We have identified a «0.15 mag. zero point error in the 
z' band, and the persistence of a low-frequency mode that 
was not removed by the flat fielding. Given these issues 
with z'-band, we cautiously analyze the {g 1 — i 1 ) CMD 
and provide photometry in Table [5] We will first analyze 
our spectroscopic results and 2MASS photometry, before 
fitting isochrones to our optical data. We will show in 
Section 14.41 that essentially identical cluster properties 
are determined from these 3 data sets. This consistency 
suggests the g' and i' zero points are likely accurate. 

The Elixir and TERAPIX pipelines were developed 
and have been successfully used to reduce similar Mega- 
Cam imaging for the CFHT Legacy Survey 24 . Despite 
these quality assurances, we consider this photometry to 
be preliminary, and we are currently working to further 
validate the accuracy of the Elixir and TERAPIX cali- 
bration. 

3. IDENTIFYING THE RUPRECHT 147 MEMBERSHIP 

We identify stars as R147 members based on their com- 
mon space motion and placement on a color - magni- 
tude diagram. Our initial membership list is drawn from 
the NOMAD and UCAC-3 astrometric catalogs and sub- 
jected to radial velocity vetting. We queried NOMAD 
and UCAC-3 for stars within a radius of 2° of the cluster 
center. Stars were accepted as candidates if their proper 
motions f ^2.1j) were within 8 mas/yr of the cluster mean 
(see Figure [3] for a proper motion vector point diag ram). 
We adopted the values of iKharchenko et al.l (|2005l) for 
the cluster center and mean proper motion: (a [h:m:s], S 
[d:m:s]) = (19:16:40, -16:17:59) and (/x a ,/x 5 ) [mas/yr] = 
(-0.6, -27.7) - after we identified the highest confidence 
members, we recalculated these locations and note no 
significant change (Table 0J. 

Table 6 gathers data for 108 stars of interest. The 
first column provides a designation internal to this pa- 
per: CWW # (CWW = Curtis, Wolfgang and Wright). 
The stars are ordered according to increasing V mag- 
nitude (provided by NOMAD). The table also includes 
the 2MASS ID (and therefore RA and Dec position); 



For a list of publications. see 

vw.cfht.hawaii.edu/Science/CFHTLS/cfhtlspublications.html 
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Figure 7. Photometric error estimates for CFHT/MegaCam 
g'r'i'z' photometry: We imaged the R147 field with four separate 
but partially overlapping pointings (see Figure [6j . We matched 
all stars in the overlap regions in our bright source catalog (short 
exposures), and find a total of 1575 unique sources with g' < 18. 
This figure plots the mean versus standard deviation of g' for the 2 
- 4 independent measurements, depending on the number of over- 
lapping regions containing the source. We find a typical value of 
A„/ = 0.035, but this is probably larger than what should be as- 
sumed for the photometric precision across the field, because one 
of the sources usually lies very close to the edge of a field, where 
the mosaic dithering is incomplete, and so the photometry is less 
reliable. 
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PPMXL proper motions in mas/yr; the preliminary 
CFHT/MegaCam g' and g' — i' optical photometry; 
2MASS J and J - K s NIR photometry; and radial ve- 
locities RVlp and RVh (Lick/Palomar and Hectochcllc, 
respectively). A membership probability is assigned to 
each of these data, according to criteria discussed below 
and summarized in Table [2] and is listed in the order: (1) 
radial distance in proper motion space from the cluster 
mean, (2) RVlp, (3) RVh, (4) proximity to cluster locus 
on the 2MASS (J~K S ) CMD and (5) CFHT/MegaCam 
{g' - i') CMD. 

The derivation of quantitative membership probabili- 
ties is precluded by the large uncertainties in proper mo- 
tion and our Lick/Palomar velocities, combined with the 
intrinsic spread in the R147 main sequence due to un- 
resolved stellar multiplicity and the possibility of differ- 
ential reddening, along with non-negligible photometric 
error. 

Instead we designate three confidence levels: 'Y' for yes 
this is consistent with cluster membership, 'P' for possible 
/ probable member, and 'N' for not likely / non-member. 
Each membership criterion is independently assessed and 
assigned a confidence level designation (whenever data 
are unavailable, a '-' is assigned instead). The following 
sections address each criterion, and establish the ranges 
for each level (summarized in Table [5]). The results from 
all fields are then reviewed and an overall membership 
confidence level is assigned to each star according to the 
same 'Y', 'P', 'N' scheme. 

We find 81 stars of highest confidence, 21 stars with 
'P' possible member status, and 6 stars with little to no 
probability of membership - at least as single star mem- 
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bers (multiple star systems could show RVs and photom- 
etry inconsistent with membership as we have defined 
it, while still being gravitationally bound members of 
R147). 

3.1. Proper motion 

We now primarily use PPMXL proper motions to as- 
sess membership, although our original membership list 
was derived from NOMAD and UCAC-3. The typical 
PPMXL errors range from 1 to 5 mas/yr in proper mo- 
tion, and so we designate stars with proper motion within 
this 5 mas/yr of the cluster mean, as 'Y' members. Only 
stars within ^8 mas/yr were considered in our initial 
candidate list, and we found that the majority of bright 
candidates had velocities consistent with cluster mem- 
bership. There are six stars that we had originally clas- 
sified as 'Y' or 'P' according to their NOMAD proper 
motions, but which have PPMXL proper motions more 
than 8 mas/yr different than the cluster's mean motion. 
Despite this large discrepancy from the PPMXL data, we 
classify these stars with conflicting proper motion data 
as 'P', and will consider their velocities and photometry 
when assigning their final probability. 

3.2. Radial velocity 

The Gener al Catalogue of Stellar Radial Velocities 
( Wilsonl [l953l) contains a single entry for a cluster mem- 
ber: HD180015 (HIP 94635, classified as KOIII, CWW 
1). iWilsonl reported RV = 41 km s~ , with quality des- 
ignation 'O, corresponding to a typical uncertainty = 
2.5 km s~ x and maximum uncertainty = 5 km s -1 . This 
was th e first and only available RV until iPakhomov et al.l 
(2009) observed three other cluster red giants: HD 
179691 (CWW 9) at 46.7 km s" 1 , HD 180112 (CWW 
10) at 40.1 km s~\ and HD 180795 (CWW 7) at 40.8 
km s , with S/N > 100, and precision estimated at 0.5 
- 0.8 km s" 1 . The RV for HD 179691 is 6 km s _1 larger 
than the other two stars, too large to be explained by 
the cluster velocity dispersion, which implies this star 
is either a SB1 binary or a non-member. We observed 
these stars at Lick/Palomar and measure RVlp = 42.1, 
41.4, 42.4 km s _1 , respectively. While our measurements 
for the second two st ars are in basic agreement with 
IPakhomov et al.l (|2009l) . the velocity for HD 179691 now 
appears consistent with the cluster mean, supporting its 
membership and corro borating its SB1 status. 

IChubak et al.l (|2012D have also measured RVs with rms 
errors ~50 m s" 1 for our 5 single-lined Keck/HIRES 
spectra (Table [3]). We list these here in km s~ , with our 
Lick/Palomar velocities in parenthesis for comparison: 
CWW 44: 41.41 (42.2), CWW 91: 41.50 (42.8), CWW 
21: 40.35 (41.9), CWW 78: 41.02 (40.8), and CWW 22: 
46.63 (51.9). We followed up CWW 22 with Hectochelle 
and measure RV = 38 km s^ 1 on both nights, and clas- 
sify it a SB1. 

Selecting the six stars above showing no evidence 
of binarity, we find a typical cluster radial velocity of 
40.86 ± 0. 56 km s -1 (if w e take only the 4 Keck stars an- 
alyzed by IChubak et all then we find the cluster RV is 
41.07 ±0.52 km s" 1 ). We take stars with RVs consistent 
with this value as high-probability cluster members. Fig- 
ure |8]plots 98 stars with Lick/Palomar velocities (shown 
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Figure 8. R147 radial velocity distribution. The gray hash desig- 
nates RVs measured from Lick/Palomar (98 stars), with RVlp = 
44 ± 3 km s — 1 . The black line plots 45 stars with Hectochelle ve- 
locities, RVh = 41.6 ± 1.5 km s — 1 . Tick marks on top indicate 
median values, color-coded to RV source: blue shows the cluster 
RV = 41.07 ± 0.52 km s _1 from the 4 Keck stars (described in 
EI3.2I1 . gray and black represent Lick/Palomar and Hectochelle re- 
spectively. The width of each distribution is as expected from the 
RV precision of each survey, and should not be interpreted as a 
resolved cluster velocity dispersion. According to the quoted RV 
precision of Chubak et al. (2012), the 0.5 km s _1 dispersion in our 
4 Keck velocities might actually be the intrinsic cluster velocity 
dispersion. 



in gray hash) with RVlp = 43.8 ± 3.2 km s . Also 
shown are 45 stars with Hectochelle velocities (black line) 
with RV H = 41.6 ± 1.5 km s -1 . The blue tick mark at 
the top shows the typical cluster velocity from above at 
40.57 km s _1 . The width of each distribution is consis- 
tent with the RV precision of each survey, and should not 
be interpreted as a resolved cluster velocity dispersion 25 
We have RVs from Lick/Palomar (RVlp) for nearly all 
stars listed in Table 6, except the 6 putative blue strag- 
glers and 4 double-lined spectroscopic binaries (SB2s); 
and RVs from Hectochelle (RVh) for 50 stars in the cen- 
tral square degree. The Hectochelle velocities are more 
precise (Figure [5]), so whenever available, RVh is used 
to determine the confidence in membership. Some stars 
have RVlp within the highest confidence interval, and 
RVh in a lower level. In these cases, if the star has Y' 
confidence level proper motions, RVlp and photometry, 
we set the overall confidence to 'P', and consider the star 
a candidate SB1 (e.g. CWW 92 has RV LP = 45 km s^ 1 
and RV H = 25 km s _1 ). 

25 M67, a much richer cluster, has a measured velocity dispersion 
of 0.5 km s from radial velocities measured by Mathicu (1983). 
Assuming virial equilibrium, the cluster velocity dispersion can be 

approximated as o\,(kms -1 ) ~ y With ~500 known mem- 

bers and similar size, we expect the M67 velocity dispersion to be 
about twice that of R147. If the 50 ms -1 RV precision estimate by 
Chubak et al. (2012) is valid, then the 0.5 km s — 1 RV dispersion 
in our Keck RVs might actually be the intrinsic velocity dispersion 
for Ruprecht 147, on par with M67 despite the lower number of 
members. 
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3.3. The Color - Magnitude Diagram and Stellar 
Populations 

Before assiging membership confidence designations, 
we check that the stars are confined to the region of color 
- magnitude space expected for a coeval stellar popu- 
lation with the properties we determine best describe 
R147. 

We mapped out this locus by simulating a rich clus- 
ter with the properties we find for R147 from isochrone 
fitting f q4.4j) . Figure [HI shows CMDs for such simulated 
clusters. In this case, we simulated a cluster with 10 6 
stars, with masses uniformly distributed between 0.6 and 
1.6 Mq. We set the binary fraction to 50%, with compan- 
ion masses uniformly distributed between zero and the 
primary mass. Differential extinction is introduced ac- 
cording to a Gaussian with /i = 0.25 and a = 0.05 mag., 
and photometric precision is set at 0.02 mag. for g'r'i'z' 
and 0.025 mag. for JHK$- The simulated photometry 
is drawn from a Padova isochrone with logi = 9.4 (2.5 
Gyr), [M/H] = +0.065, and m — M = 7.35 mag. We 
bin, log-scale, and smooth the synthetic photometry to 
highlight the R147 locus in color - magnitude space. 

Stars overlying the shaded region (basically, the re- 
gion bound by the single star and equal mass binary se- 
quences) are given the highest confidence designation. 
The simulation demonstrates that atypical differential 
reddening along a particular line of sight or relatively 
high photometric error can place stars outside the locus. 
Stars in these regions are assigned 'P'. These could also 
be triple systems or exotic products of stellar mergers. 
CWW 67 is the only star existing beyond the equal mass 
triple sequence, and we assign it the lowest designation, 
'N'. 

Confidence assignment is an iterative process, since 
we identify high-confidence members using isochrone fits, 
and these fits require a list of high-confidence members 
so that unlikely or non-members do not throw off the fit. 

3.3.1. Stellar Populations 

Blue Stragglers: In addition to the potential triple 
systems, 5-6 stars occupy a space of the CMD out- 
side the cluster locus beyond the main sequence turnoff 
(MSTO): 6 clearly separate in the 2MASS CMD, but 
only 5 in g'r'i'z' . These five stars have proper motions 
consistent with the cluster, but lack RV measurements 
due to rotational line broadening (CWW 24, the sixth 
outlier in 2MASS, does have a measured RVlp = 41.8 
km s , and so we assume that 2MASS photometric er- 
ror is responsible for scattering it out of the cluster lo- 
cus). We classify these 5 stars as blue stragglers (see 
Table 6, blue stragglers are listed as 'BS' in the Notes 
column). For the photometric probabilities, instead of 
the 'Y/P/N' scheme, we assign a 'B' for blue straggler. 

Red Giants: We find 11 red giants in the cluster. The 
TERAPIX photometric errors suggest that only the four 
brightest red giants are saturated in g'r'i'z' even with 1 
second exposures. Other stars down to «9.5 are quoted 
as saturated in each band across the 4 fields. After con- 
sulting the raw frames, reduced images, and considering 
the 9.5 magnitude saturation limit found for other stars, 
we conclude that the entire red giant branch has unre- 
liable optical photometry. This explains the apparent 
mismatch between our best isochrone fit and the optical 



RGB. 

CWW 14 is fainter and has reliable photometry, al- 
though it is 0.15 mag. blueward of the red giant branch 
in (q' — i'), an d 0.08 mag. (4-cr) blueward in (J — Ks). 
iMathieu etafl (|1990f ) identify a SB1 system in M67, 
S1040, which lies 0.2 mag. to the blue of the red giant 
branch in (B — V) . This system was previously suggested 
to consist of a star further down the giant branch with a 
companion star near the MSTO. [Landsman et al.l (1997) 
identified broad Lyman absorption features, demonstrat- 
ing that the companion is actually a hot white dwarf, and 
that the system likely underwent a period of mass trans- 
fer. CWW 14 is an outlier in NIR, making it less likely 
to be a "red straggler" , and is probably a MSTO - RGB 
binary. 

Main Sequence dwarfs: We use our best isochrone 
fit to determine approximate spectral types for the R147 
membership. We assu me masses of 1. 1 M for GO and 
0.8 M for K0 dwarfs (jZombeckl 120071) . then locate the 
boundaries in the CMD from the isochrone. We find that 
the MSTO is located around mid-F. The subgiant branch 
down to F8 on the main sequence is well populated with 
«52 stars. We also identify «27 G dwarfs and «8 K 
dwarfs down to mid-K (We quote approximate numbers 
because of the approximate nature of our spectral typ- 
ing). The 9 stars lacking g'r'i'z' photometry appear in 
the 2MASS CMD as follows: 2 red giants, 3 MSTO F 
stars, 2 G dwarfs and 2 K dwarfs. 

This method ignores the existence of binaries but illus- 
trates the top-heavy nature of our membership list. This 
is likely due to a combination of observational bias and 
cluster evaporation (star clusters tend to lose their low- 
est mass members first, and 'evaporate' from the bottom 
up). The typical NOMAD proper motion error is ~10 
mas/yr by V ~ 12. The K dwarfs have V > 13, mak- 
ing candidate identification from proper motions diffi- 
cult. Therefore we are almost certainly missing signifi- 
cant numbers of low mass dwarfs. 

3.4. Notes on Particular Stars 
3.4.1. Apparent Non-members 

CWW 77 has RVs inconsisent with R147 (RV L p = 51 
km s~\ RV H = 52 km s _1 ), but its (g'-i') and (J-K s ) 
photometry place it near the equal mass binary sequence, 
so it could be an SB1. 

CWW 67 has a low RV L p = 34 km s^ 1 and is 1.3 mag. 
above the (g' — i') main sequence and 1.4 mag. above in 
(J — Ks), but an equal mass triple would sit 1.2 mag. 
above the main sequence, so membership seems very un- 
likely (although a fourth lower mass companion could 
theoretically explain these discrepancies, so membership 
is difficult to definitively rule out). 

CWW 72 sits wl mag. above the main sequence in 
(g' — i') and w0.75 mag. in (J — Ks) - if a member, 
this could be a triple or an equal mass binary with inac- 
curate optical photometry. CWW 72 is a SB2, as seen 
in the Keck/HIRES spectrum. This star was observed 
previously at Lick and the spectrum exhibited no sign of 
binarity (otherwise we would not have selected it for ob- 
servation with Keck) . We also observed CWW 72 on two 
consecutive nights with MMT/Hectochelle. The CCF 
from the first night exhibits a tall and sharp peak with 
RVhi = 46.6 km s _1 . The CCF from the second night is 
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Figure 9. Each panel illustrates a source of main sequence broadening, demonstrating why the R147 main sequence might appear thicker 
than a textbook "beads on a wire" CMD. The R147 stars are plotted in green and saturated stars are plotted with an open circle. A Padova 
isochrone is overlaid in red with age = 2.51 Gyr, m — M = 7.32 (d = 291 pc), Ay = 0.23, and [M/H] = +0.065. The shaded regions in 
each panel represent simluations of 10 6 stars, with masses uniformly distributed between 0.06 and 1.6 Mq, and photometry queried from 
the previously quoted Padova model. The simulated photmetry has been binned (0.005 mag. in g' — i' , 0.01 mag. in g'), log scaled, and 
smoothed with a 5 pixel boxcar, to highlight possible regions of color-magnitude space occupied by R147 members. The top-left panel only 
includes photometric error, set at cr g r ir = 0.02 and assuming normally distributed errors. The top-right panel only includes binaries, with 
the binary fraction set at 50%, and the secondary masses uniformly distributed between zero and the primary mass. The bottom-left panel 
only includes differential extinction, normally distributed about the typical cluster value of Ay = 0.25, with SAy = ./V(0,0.05). 



lopsided, suggesting that the signature of the companion 
was beginning to manifest and that the period of this 
system could be on the order of days. The CCF shape 
and resulting radial velocity from the first night point to 
a systematic velocity «5 km s _1 greater than the R147 
bulk motion, which cannot be explained by the cluster's 
velocity dispersion or RV precision. If CWW 72 is a 
member, then it is (at least) a triple, perhaps with two 
approximately equal mass primary components orbiting 
with a period of days, and a fainter companion modulat- 
ing the RV on a longer timescale (needed to explain the 
5 km s _1 systematic offset). 
Finally, CWW 50 sits 0.05 mag. blueward of the 



(g' — i') main sequence, but is on the (J — Kg) main 
sequence. In Section 14.11 we discuss the possibility that 
this star is less extinguished and reddened than the rest 
of the cluster. If this is not the case, perhaps a hot white 
dwarf is pulling it blueward while not introducing much 
NIR flux, or there is an atypically large photometric er- 
ror in one of the optical bands (2-cr), or else CWW 50 is 
not a member. We list it as 'P' because it is only incon- 
sistent in (g' — i'), and while RVlp = 47.6 km s _1 , the 
Lick/Palomar velocity precision does not rule out mem- 
bership. 

3.4.2. Notes on 2MASS photometry 
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Table 2 

Criteria for membership 



Data Source Highest: Y Probable: P Low or Non-member: N 

r M a NOMAD, UCAC-3, Adam Kraus < 5 5-8 > 8 

Radial Velocity 13 Lick & Palomar 39 -47 36 - 39, 47 - 50 33 - 36, 50 - 53 

Radial Velocity Hectochelle 40 - 43 38.5 - 40, 43 - 44.5 else 

(J — Kg) CMD 2MASS Overlaps with simulation ± 0.2 mag. beyond equal mass triples 

(g 1 — i') CMD CFHT/MegaCam Overlaps with simulation ± 0.2 mag. beyond equal mass triples 



-> 10 



Note. — Sec Section^ for a discussion of membership criteria. Values in parenthesis denote acceptable ranges. Values equal to endpoints 
are assigned to the higher level 

a Radial distance in proper motion space from the mean value for R147, with units of mas/yr. 

Radial velocities measured in km s — 1 . 
c Sec Section 13 .3! and /or Figurc[9]for discussion 

Our g'r'i'z' imaging shows CWW 51 is an optical dou- 
ble, with a star 1.65 arcseconds away with a similar 
g'r'i'z' SED (the magnitude difference in each band is 
0.02 - 0.05 mag. between the two stars). The com- 
ponents are separated by 1.65", which translates into a 
minimum physical separation of 495 AU, assuming a clus- 
ter distance of 300 pc. This suggests that the pair actu- 
ally form a wide binary, although their angular proximity 
could also be explained by a chance alignment. This dou- 
ble was not resolved in the 2MASS Point Source Catalog. 
Adding 0.75 mag. to the J band magnitude (halving the 
brightness, to reflect just the one star) moves CWW 51 
next to the stars it neighbors in the (<?' — i') CMD. Table 
6 quotes, and the figures in this work plot, the 2MASS 
photometry for CWW 51, despite this realization, al- 
though we do include a footnote referencing this in the 
Table. 

The 2MASS Point Source Catalog provides PSF pho- 
tometry by default in most cases. Figure [10] shows 23 
outlier stars on either side of the {J —K$) main sequence, 
out of the 80 stars 'Y/P' stars with aperture photome- 
try that are not blue stragglers. If aperture photometry 
is used instead, each of these stars moves towards the 
cluster locus. No star already in the locus shifts appre- 
ciably outside when aperture photometry is used instead. 
One star slides from blueward from the main sequence 
by 0.037 mag., or approximately equal to the J and Ks 
errors added in quadrature. The fact that the majority 
of outliers' photometry systematically moves toward the 
cluster locus suggests to us that for many stars in these 
fields and at these magnitudes, the aperture photometry 
is superior. We do not assign lower confidence levels to 
PSF photometry outliers, if their aperture photometry 
is consistent with membership. We include the aperture 
photometry for 18 stars in Table 6 and all other figures 
(CWW 22, 24, 27, 28, 30, 37, 38, 43, 48, 49, 57, 59, 67, 
68, 90, 91, 94, and 100). 




Figure 10. PSF photometry (green circles) vs aperture photom- 
etry (red triangles) for 18 outliers on the 2MASS (J - K s ) CMD. 
All shift closer to the locus when the aperture photometry is used 
instead of the default PSF photometry. No star already in the lo- 
cus shifts appreciably outside when aperture photometry is used 
instead. One star slides from blueward from the main sequence by 
0.037 mag., or approximately equal to the J and Kg errors added 
in quadrature (not shown). 



mass status and the validity of our isochrone fit. 
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3.4.3. SB2 systems 

CWW 64, 65, 66, 68, and 72 showed double-peaked 
cross-correlation functions in one of the RV epochs, in- 
dicating these systems to be nearly equal mass binaries 
(the case of CWW 72 is discussed above). FigurelTTIplots 
the (g' - i') and (J - K s ) CMDs, with the five SB2s 
highlighted red. The single star and equal mass binary 
sequences from our best isochrone fit (Padova model) are 
plotted in green. All five SB2s are clustered around the 
equal mass binary sequence, corroborating their equal 



3.4.4. SB1 binary candidates from discrepant RVs 

Figure [TT1 also plots in cyan 10 stars with Hectochelle 
RVs inconsistent with the cluster: CWW 19, 22, 27, 53, 
69, 70, 77, 92, 95, 99, and 106 (CWW 19 is a red gi- 
ant. We do not include it in the optical plot because 
it is saturated). In Section 13.21 we suggested CWW 
92 is a SB1 from the 20 km s _1 difference between RV 
epochs. All other SB1 candidates have similar RVs from 
Lick/Palomar and Hectochelle, and are considered can- 
didate SBls because RVh is at least a few standard de- 
viations away from the cluster average, although this 
discrepancy forces a 'P' classification (except CWW 77, 
which we classify as 'N' as noted above). 

26 CWW 68 actually sits on the equal mass triple sequence. Both 
CWW 64 and 65 are midway between the equal mass binary (-0.75 
mag.) and triple (-1.2 mag.) sequences, at 1.0 mag. brighter than 
the main sequence, which can occur when an equal mass binary 
system, which manifests as the SB2, has a third companion with 
50% the luminosity of each primary. 
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Figure 11. Ten stars with discrepant Hcctochelle RVs are plotted in cyan (RV more than 1.5 km s away from cluster average). Five 
SB2s are plotted in red. The rest of the cluster is shown in gray, with our best fit isochrone (Padova model) overplotted in dark blue 
showing the single and binary sequences (logt = 9.4, [M/H] = +0.064, m — M = 7.32, Ay = 0.23). All of the SB2s were identified by 
a double-peaked cross-correlation function, which indicates these should be nearly equal mass ratio systems. The fact that they are all 
clustered around the binary sequence (shifted "up" 0.75 mag.) coroborates their equal mass status and the validity of our isochrone fit. 
All but 3 SB1 candidates are also shown near or on the binary sequence. The square shows CWW 99 on the (</ — i') single star sequence, 
but in (J — Kg) the star is 0.5 mag above this sequence. The other two stars might be high-mass ratio systems (and therefore do not 
manifest in shifts on the CMD), stars that have received gravitational kicks (so their RVs are no longer consistent with the cluster), or are 
non-members. 



Figure HT1 shows all but CWW 70 are clustered around 
the binary sequence. CWW 70 has RVlp = 38 km s" 1 
and RVh = 39 km s _1 . Despite this 2-er discrepancy, we 
cannot rule out membership. The star might have a low 
mass companion or received a gravitational "kick" . For 
example, a G2 V (1 M ) - M0 V (0.5 M ) binary system 
with semi-amplitude velocity K\ = 5 km s~ , zero ec- 
centricity and zero inclination, will have a period P « 12 
years. The large luminosity difference in a high mass ra- 
tio binary means the system will not stand out in either 
photometry (it will lie on the single star sequence) or in 
spectra (the secondary is too faint to manifest as a SB2). 
This means although the RVs for CWW 70 are incon- 
sistent with single star membership, we cannot rule out 
the possibility of a low mass companion at large separa- 
tion, which would induce a measureable velocity offset, 
but modulated at a period much longer than our 2 year 
baseline. 

CWW 99 is plotted as a blue square in Figure QTJ It 
sits on the single star (<?' — i') sequence, but is 0.5 mag. 
above the (J — Ks) sequence. This can be explained 
by a low mass companion, which would show up more 
prominently in NIR than optical. 

4. INFERRING CLUSTER PROPERTIES 

The properties of stellar clusters (age, composition, 
distance and interstellar extinction) are commonly es- 
timated by fitting isochrones to broadband photomet- 
ric color - magnitude diagrams (CMD). Often a "chi- 
by-eye" technique is employed, where sets of isochrones 
representing varying cluster parameters are overlaid on 
a CMD, and the apparent best fit or series of best fits 
are selected to establish acceptable ranges for these fun- 
damental parameters. This technique can be successful 



when one or more of these properties can be well con- 
strained. For example, clusters might be nearby or sit 
above/below the galactic plane and suffer little extinc- 
tion, and the closest benchmarks have parallax distances 
from Hipparcos and/or HST/FGS. These independent 
constraints break the high degree of degeneracy between 
each variable (e.g. metallicity works in a similar direction 
to the interstellar reddening vector). 

None of Ruprecht 147's properties have been previ- 
ously well measured. At a distance of over 200 pc, the 
cluster has a HIP2 distance measurement from 3 stars 
that appears unreliable and is apparently too close by a 
significant fraction ( fl4.7[) . We first describe our efforts to 
ind ependently constrain the interstellar extinction with 
the lSchleeel et all |l998) dust map ( fc|4.1[) , and the com- 
position from spectroscopic analysis f £|4.2[) , then we will 
use isochrone models to determine the cluster's age and 
distance. Specifically, we fit a spectroscopically derived 
T e ff - log g diagram with Padova isochrones with abun- 
dances fixed by the spectroscopic metallicity ( fc|4.3[) . We 
then query the resulting best fit Padova isochrone for 
a star with T e s and log g closest to the values for the 
early G dwarf we derive with SME, and perform a brute 
force x 2 SED fit for distance and visual extinction to the 
resulting synthetic g'r'i 1 JHK$ photometry. 

Next we develop an efficient and automated 2D cross- 
correlation isochrone fitting techinque, and fit Padova 
and Dartmouth models to the (g' — i') and (J — Ks) 
CMDs to determine the age, distance and visual extinc- 
tion. We will see that the NIR and optical fits agree very 
well with each other (Figures [T71 [T5| . We also fit the 
(g 1 — i') CMD with the r 2 maximum-likelihood method, 
and derive consistent results, validating our isochrone 
fitting technique. We find an age from both NIR and 
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Figure 12. Visual extinction map with R147 members de signated 
by red dots. Ay calculated from the dust map of ISchlegel et afl 
(1998), assuming Ry = 3.1. The map was smoothed with a 3 
arcminute boxcar. 



optical photometric isochrone fitting consistent with our 
spectroscopic results, which supports our earlier decision 
to break the degeneracy between these parameters with 
this age (03J). 

We will compare results from Padova, Dartmouth and 
PARSEC models; and then discuss the differences be- 
tween our photometric distance and the parallax dis- 
tance inferred from 3 H1P2 cluster members. Finally, 
we synthesize the results from these various isochrone 
fits and present our preferred set of parameters describ- 
ing the age, composition, distance and visual extinction 
for Ruprecht 147 (©. 

4.1. Interstellar Extinction 

The large apparent size of R147 on the sky introduces 
the possibility of differential extinction acr oss the cluster. 
Figure [T2l plots Ay from the dust map of ISchlegel et al.l 
(1998), assuming Ry = 3.1 27 . Many cloud structures are 
apparent in the field, showing Ay to vary from 0.3 - 0.5 
mag. 

Considering the cluster's proximit y, some of this dust 
undou btedly lies beyond the cluster. iDrimmel fc Spergel 
(2001) determine the Sun to lie 14.6 ± 2.3 pc above the 
Galactic midplane, and measure a dust scale height of 
hd = 188 pc at the Solar Circle. The Galactic latitude of 
R147 ranges from -12° to -14°. At a distance of 250 - 300 
pc, this latitude places R147 50 - 70 pc below the Sun, 
and 35 - 55 pc below the midplane (less, if a larger Z Q 
is assumed), or about 20-30% of the dust scale height. 
The local bubble has very little dust in it out to ~ 150 

27 Many studies have demonstrated that on average, Ry = 3.1 
for a diffuse ISM. This relationship between extinction and redden- 
ing is of course dependent on the composition and physical condi- 
tions in the intervening ISM. Considering Ruprecht 147, there are 
no dense molecular clouds along the line of sight. We also do not 
expect any additional reddening intrinsic to the cluster or its stars, 
because this is an old cluster, and it and its stars are no longer en- 
shrouded in dust and gas. Recently, IJones et al. ( 2011) used over 
56,000 Sloan M dwarf spectra to map visual extinction and Ry in 
the nearby Galaxy. While their Ry distribution peaks near 3.1, 
their median value is actually 3.38, and ranges from 2 to 5.5. We 
will assume the canonical value of Ry = 3.1 for the diffuse ISM 
for this preliminary study and intend to return to this issue in a 
future work. 



pc (|Lallement et al.|[2003^ . suggesting it is possible that 
most of t he d ust is behind the cluster. 

Section 13.31 described star cluster simulations used to 
assess photometric membership probabilities, by intro- 
ducing photometric scattering sources to explain the ob- 
served width of the R147 main sequence, including pho- 
tometric error, binarity, and differential extinction. Fig- 
ure [9] plots simulated cluster CMDs including each of 
these photometric scattering sources separately, and all 
combined. The binarity simulation demonstrates that 
the single star and equal mass binary sequences pinch to- 
gether near the MSTO. Differential reddening smears the 
CMD along a negative-sloped diagonal (extinction plus 
reddening) . If there is non- negligible differential redden- 
ing, the "binary pinch" should be smeared out. 

Unfortunately, we have only identified 6 members at 
this pinch. Figure [13] shows four are confined within the 
pinch. CWW 50 sits 0.08 mag. blueward, and CWW 
53 sits 0.06 mag. redward. These values are 2-3 times 
larger than the expected photometric error. CWW 53 
has RVs from two epochs at 5 km s™ 1 greater than the 
cluster mean, and photometry placing it near the equal 
mass triple sequence. The uncertainty in membership 
and multiplicity means we cannot use CWW 53 to test 
for differential reddening. CWW 50 must be reddened 
by 8 Ay = 0.13 mag. in order to place it on the single 
star sequence. In Section 14.31 we find an optimal Ay = 
0.23 mag., so a particular line of sight of Ay = 0.10 
mag. is not impossible. Ideally, we would like to check if 
the nearest R147 neighbors to CWW 50 also appear less 
extinguished than expected. Unfortunately, the nearest 
neighbor is ~ 1 0' aw ay. It is also noteworthy that the 
ISchlegel et al.l (|1998[ ) dust map extinction along this line 
of sight is Ay = 0.296 mag., the lowest value in the entire 
field, and 0.13 mag. lower than the median Ay for the 
region of radius r =1°.2 encompassing all R147 members. 

The R147 main sequence is thicker than 1 magnitude 
at various points. Unfortunately in many cases, our ra- 
dial velocities do not have sufficient precision to firmly 
establish these stars as cluster members (we still desig- 
nate them 'Y' members because the RVs are consistent 
with the R147 bulk motion, within the precision of our 
Lick/Palomar RV survey). We also only have one RV 
epoch for the majority of stars, and so stellar multiplic- 
ity is impossible to diagnose at this point. More precise 
velocities are required before we attribute these photo- 
metric outliers to differential extinction. 

There is also a strip of 7 main sequence dwarf stars 
blueward of our best isochrone fit in the (g' — i 1 ) CMD. 
Adding 8 Ay — 0.05 mag. to these stars shifts them onto 
the isochrone, This translates into a 0.03 mag. shift in 
color, which is within the photometric precision, and so 
differential extinction is not required to explain these 
stars apparent blueward offset in the optical CMD. No 
net offset is seen in the NIR CMD. 

We will postpone further investigation into differential 
extinction to a future study, and in this work will fit 
single Ay models, with values constrained by the dust 
map at Ay < 0.5. 

4.2. Metallicity 

We analyzed five Keck/ HIRES spectra with Spec - 
troscopy Made Easy (SME, iValenti fc Piskunovl Tl 996). 
using the procedure described in IValenti fc Fischen 
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Figure 13. Looking for differential extinction. The 'Y' and 'P' 
members are plotted with 0.03 mag. error bars in color and magni- 
tude. Single star (yellow) and equal mass binary (cyan) sequence 
isochrones (Padova) are plotted with logt = 9.4, [M/H] = +0.1, 
m — M = 7.48, and Ay = 0.23. The gray shading illustrates the 
results of our Monte Carlo cluster simulation, including binaries 
and neglecting differential extinction and photometric error, from 
Figurc[9](top - right panel). The single and binary sequences cross 
at g' 10.9, forming a "pinch" at the turnoff. Differential extinc- 
tion should smear this pinch along the reddening vector shown in 
blue. The three outliers redward of the binary sequence could be 
triples. CWW 50, the one outlier blueward of the single star main 
sequence, might either suffer atypically large photometric error, or 
is not a member, or else needs to be reddened by SAy = 0.13 in 
order to place it o n the main sequence . Interestingly, the visual 
extinction from the Schlcgcl ct al. ( 1998) dust map at the position 
of CWW 50 is the lowest in the field, and exactly 0.13 mag. lower 
than the median value for the region of radius r =1°.2 centered on 
the cluster and encompassing all R147 members. 

(2005). SME uses the Levenberg-Marquardt algorithm 
to fit observed echelle spectra with synthetic spectra 
generated assuming LTE and plane-parallel geometry, 
yielding effective temperature, surface gravity, metallic- 
ity, projected rotational velocity, and abundances of the 
elements Na, Si, Ti, Fe, and Ni 28 

Upon obtaining an initial fit to a spectrum, T e g- was 
perturbed ±100K and run again. The three solutions 
were then averaged, with the standard deviation set as 
the parameter uncertainty, except in cases where this un- 
certaint y is less than the s tatist ical uncertainities mea- 
sured in lValenti fe Fischerl (|2005l ): 44 K in effective tem- 
perature, 0.03 dex in metallicity, 0.06 dex in the loga- 
rithm of gravity, and 0.5 km s _1 in projected rotational 
velocity. Additional correction s are applied to the fina l 
values based on the analysis in lValenti fe Fischerl (2005) 
of Vesta and abundance trends in binary pairs with T e g 
(see Figure [14] for the CWW 44 Keck/HIRES spectrum 

28 Valenti k. Fischer (2005) did not solve for magnesium abun- 
dance because of the degeneracy between [Mg/Fe] and log g when 
fitting the synthe t ic spectra to the gravit y-sensitive Mg b triplet 
(Thackerav 1939). Fuhrmann ct al. (1997) recommends fitting for 
[Mg/Fe] first, using the weak Mg I lines at A4571 or A5711, and 
then fixing the abundan c e and fitting the Mg b wings to derive 
logg. Valenti & Fischer (2005) decided to exclude from analysis 
wavelengths A < 6000 A, except for the region encompassing the 
Mg b triplet, because of severe line blending in the blue in cool 
stars which would complicate the spectral synthesis fit. 




5165 5170 5175 5180 
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Figure 14. A Keck/HIRES spectrum of the order encompassing 
the Mg b triplet for CWW 44, an F MSTO star, is shown in black. 
The synthetic spectrum resulting from our SME analysis is overlaid 
in red. The spectrum segments included in the fit are highlighted 
in purple, and the salmon stripes identify the continuum regions. 
Table [3] lists our SME results for this and four other stars. 



and SME synthetic spectrum fit in the order encompass- 
ing Mg b). 

Our results for the five stars are presented in Table [3] 
and indicate that the cluster has a slightly super- Solar 
metallicity of [M/H] = +0.07 ± 0.03, from CWW 91, 
44 and 21 29 . We neglected the results from CWW 22 
because it is hottest (complicating the fit to gravity, de- 
scri bed in the next su b sectio n) and has the poorest xl 
fit. IValenti fe Fischerl (|2005|) suggest using [Si/H] as a 
proxy for the a-process abundance. We find [a/Fe] = 
[Si/H] - [Fc/H] « 0.0 ([Si/Fe] = -0.03 and 0.0 for CWW 
44 and 91). 

We find a much lower metallicity for CWW 78, [M/H] 
= -0.11 ± 0.03 and [Fe/H] = 0.0 ± 0.02. This outlier has 
otherwise satisfied every criterion for membership, with 
proper motions, photometry and a precise RV all consis- 
tent with the cluster. For this work, we will assume that 
this peculiar metallicity can be explained by a compli- 
cation in the SME analysis, and will look into this in a 
f uture study. 

IPakhomov et al.l ([2009) analyzed high-resolution, high 
signal-to-noise spectra of three red giant members of 
R147, and their results are compiled in Table [3] The 
first star, HD 179691, has a radial velocity inconsistent 
with the cluster, indicating it is either an SB1 or not a 
member. The other two red giants show super-Solar iron 
abundance, consistent with our SME results. 

In the future, we will more rigorously determine the 
cluster metallicity combining photometry, spectroscopy, 
and cluster properties. 

4.3. Fitting Isochrone Models to Spectroscopic 
Properties 

Figure IT5l shows the T c g - logg diagram resulting from 
our SME analysis, along with a (</ — i') CMD for the 5 
stars, along with their CWW IDs. The CMD shows that 
our SME results place the stars on the T e ff _ log g diagram 
with the correct relative positions. We have selected 

29 We will often use [M/H] = +0.064 throughout this work, 
because of the way the Padova isochrone webtool parameterizes 
metallicity: [M/H] = log(Z/Z0), with Z = 0.019. Setting Z = 
0.022 gives [M/H] ?s +0.064. 



Table 3 

Spectroscopic analysis of R147 stars 



Property 


CWW 44 a 


CWW 91 a 


CWW 21 a 


CWW 22 a 


CWW 78 a 


HD 179691 b 


HD 180112 b 


HD 180795 b 


Type 


mid-F MSTO 


GO/2 V 


subgiant 


mid-F MSTO / SB1? 


late-F V 


Kl III 


K0 III 


K0 III 


T cS (K) 


6273 (5) 


5747 (62) 


6129 (25) 


6350 (80) 


6115 (52) 


4573 (80) 


4733 (80) 


4658 (80) 


logg (gm cm s~ 2 ) 


4.11 (0.02) 


4.35 (0.11) 


3.79 (0.07) 


3.6 (0.06) 


4.27 (0.08) 


2.28 (0.15) 


2.53 (0.15) 


2.43 (0.15) 


RV (km s" 1 ) 


41.41 


41.50 


40.35 


46.63 


41.02 


46.7 


40.1 


40.8 


vsini (km s _1 ) 


6.87 (0.69) 


0.32 (0.33) 


6.50 (0.61) 


6.91 (0.73) 


6.09 (0.65) 








[M/H] 


0.07 (0.01) 


0.06 (0.03) 


0.09 (0.03) 


-0.01 (0.04) 


-0.11 (0.03) 








[Na/H] 


0.22 (0.02) 


0.23 (0.01) 


0.12 (0.01) 


-0.02 (0.08) 


-0.14 (0.03) 


0.24 


0.16 


0.24 


[Si/H] 


0.14 (0.01) 


0.11 (0.02) 


0.10 (0.01) 


0.00 (0.03) 


0.02 (0.02) 


0.08 ± 0.08 


0.15 ± 0.06 


0.25 ± 0.08 


[Ti/H] 


0.28 (0.02) 


0.16 (0.02) 


0.25 (0.03) 


0.03 (0.05) 


0.06 (0.04) 


-0.03 ± 0.07 


0.04 ± 0.05 


-0.02 ± 0.04 


[Fc/H] 


0.17 (0.01) 


0.11 (0.02) 


0.22 (0.01) 


0.08 (0.03) 


-0.00 (0.02) 








[Fci/H] 












0.03 ± 0.06 


0.14 ± 0.06 


0.16 ± 0.04 


[Fell/H] 












-0.02 ± 0.09 


0.07 ± 0.07 


0.08 ± 0.04 


[Ni/H] 


0.07 (0.01) 


0.05 (0.02) 


0.13 (0.01) 


-0.03 (0.06) 


-0.02 (0.03) 


-0.04 ± 0.08 


0.04 ± 0.06 


0.12 ± 0.07 


X? 


2.90 


2.86 


4.07 


5.88 


2.18 








logt (years) 












9.2 ± 0.5 


9.0 ± 0.4 


9.0 ± 0.4 



Note. — Rows and SME statistical uncertainties: (Type) Rough Spectral Typo, (T„gr) Effective Temperature: a — 44 K, (log g) Surface Gravity: er 
= 0.06 dex, (RV) Radial Velocity (v sin?) projected rotational velocity: a = 0.5 km s -1 , ([M/H]) Metallicity = log 1Q Z/Z Q ([Na/H] .. [Ni/H]) Sodium, 
Silicon, Titanium, Iron and Nickel abundance: a — 0.03 dex, (x«) Reduced x 2 or " the fit, (log t) age in years assuming d — 280 pc 
a Our SME analysis results of Kcck/HIRES spectra 

Red giants analyzed by [Pakhomov ct al. (2009), reproduced here for comparison. 
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Padova isochrone models with [a/Fe] = and [M/H] = 
+0.1, and attempt a "chi-by-eye" fit by overlaying mod- 
els with ages at ~2, 2.5, and 3~Gyr (logt = 9.3,9.4,9.5). 
Fitting the T e g - logg diagram is powerful because it 
is independent of m — M, Ay , and color - temperature 
transformations. 

The MSTO star CWW 44 in theory provides a tight 
constraint on age and metallicity, but we are cautious of 
the accuracy of logg, because the broad wings of the Mg 
b lines provide the gravity constraint, and their sensi- 
tivity decreases at higher temperature and lower gravity. 
Jeff Valenti has suggested that the Mg b wings provide 
useful constraints on gravity for dwarfs coller than ~6200 
K 30 , which is app roximately the tempera ture for 4 of 5 
stars we analyze. iFuhrmann et al.l (|1997f) is able to de- 
rive an accurate logg for Procyon (F5V) from Mg b, so 
perhaps our concern is unwarranted. Assuming we have 
derived accurate stellar properties, we find that models 
fit best with logt = 9.4 ± 0.05 (2.25 - 2.8 Gyr), which 
encompass the error bars of CWW 44. 

The models barely pass through the error bars for the 
early G dwarf, CWW 91, even though this should be 
the one star of the five we know has broad enough Mg 
b wings to provide adequate constraint on logg, since it 
is most similar to the Sun. If we fix logg according to 
the Padova isochrone (logt = 9.4 (2.5 Gyr) and [M/H] 
= +0.064), we find logg = 4.45 instead of 4.35, at mass 
M=1.O3M . 

We queried the g'r'i' and JHK$ photometry for a 1.03 
Mq star from this isochrone and perform a brute force 
least-squares fit for distance modulus and visual extinc- 
tion for CWW 91 in the range m - M = 7-8 and 
Ay = — 0.5, with 0.01 mag. step sizes, JHKs errors 
according to 2MASS {a = 0.023, 0.026, 0.021 mag.) and 
g'r'i' errors set to a — 0.03 mag. We find a minimum y 2 
at m - M = 7AA ± 0.02, Ay = 0.25 ± 0.03. (see Figures 
IT5l H6]> . We perturbed [M/H] ±0.02 dex (our error bar 
from the SME analysis) and logt ±0.5 (our error bar from 
fitting isochrones to the T e g - log g diagram) , then re- fit 
and find uncertainties of 0.04 and 0.01 mag. for m — M 
and A v . We then perturbed T cS by ±50K (the SME 
statistcal error bar) and re-fit, and find uncertainties of 
0.06 and 0.04 mag. for each parameter. 

Adopting these conservative errors, we find m — M = 
7.44 ± 0.06 and Ay = 0.25 ± 0.04. We repeated this 
analysis with solely the optical g'r'i' photometry, and 
then with just the 2MASS NIR JHK S photometry, and 
find m-M = 7.44 and 7.41 respectively, and Ay = 0.18 
for both cases. This give us confidence in the g' > 10 
CFHT photometry. 

We introduced reddening to the synthetic SED using 
the relationships provided by the Padova CMD web- 
site 31 : A g ,/A v = 1.167, A r ,/A v = 0.860, A v /A v = 
0.656, Aj/Ay = 0.290, A H /A V = 0.183 , A K JA V = 
0.118. 

4.4. Fitting isochrone models to broadband photometry 

Up to this point, we have fit isochrones to a T e s - 
logg diagram with values derived with SME for 5 stars 
with Keck/HIRES spectra. Now we perform a more tra- 

30 http:/ /www-int. stsci.edu/~valenti/sme. html 

31 for a G2 dwarf, using a Cardclli et al. (1989) extinction curve 
with Ry =3.1 



ditional fit to the broadband optical and NIR photom- 
etry of all cluster members, and find results consistent 
with our spectroscop ic solution. We in i tially worked 
with Padova models (jGirardi et al.l 120001 : iMarigo et al.1 
12008ft 32 because they were the only group, to our knowl- 
edge, that provided isochrones in the CFHT/MegaCam 
g'r'i' z' filter set. Aaron Dotter has since released Dart- 
mouth models in this set, and we will briefly compare 
results bewteen these two models in Section FOH 

Given the high degree of degeneracy between age, dis- 
tance and visual extinction (we fix metallicity according 
to our SME result, [M/H] = +0.064, [Fe/H] = +0.1), we 
developed an automated isochrone fitting technique that 
can efficiently co ver a large parame t er spa ce. Inspired by 
the r 2 method of Naylor & Jeffries (2006) , described and 
utilized in Section [4.51 our method simulates a star clus- 
ter with a particular age and composition, and computes 
2D cross-correlations between the resulting synthetic and 
actual CMD density distributions in order to find the 
distance and visual extinction that best aligns the model 
to the data. The age and composition control the mor- 
phology of the stellar locus on a CMD, while changes 
in distance modulus and extinction simply translate the 
locus across the CMD plane. 

Using our cluster simuator ( fl3.3[) . we map the stellar 
locus with 10 4 stars, the binary fraction set at 70%, and 
no differential extinction, for ages ranging from 1 to 4 
Gyr. The simulated cluster CMD is binned by 0.005 
magnitudes in each dimension, as is the actual R147 
CMD. Photometric error is not introduced to the syn- 
thetic photometry. Instead, the position for each star on 
the real CMD is broadened by a 2D gaussian according 
to the assumed photometric error for each band: we use 
0.02 mag. for g' and i' , and the errors provided by the 
2MASS Point Source Catalog for J and K s . 

For the NIR (J — Ks) fit, all 'Y' members were binned 
except for the blue stragglers, and the giants that appear 
to be undergoing core helium fusion according to their 
optical and NIR photomet r y offs et from the main red 
giant branch. IDotter et al.l (|2008[ ) note that the Padova 
isochrones are hotter in the lower main sequence, starting 
at w 0.8 M Q , than other commonly used models includ- 
ing the DSEP (Dartmouth Stellar Evolution Program) 
and Yale-Yonsei (Y 2 ). The left and central panels of 
Figure [501 illustrates this difference in T e g - log g and 
(g' — i'), but the right panel shows that the discrepancy 
does not significantly affect the NIR isochrone. We do 
remove the K dwarfs from our optical isochrone fits. We 
also discard the optical red giant branch because of the 
saturated i' band photometry. 

We compute the 2D cross-correlation between the real 
and synthetic CMD distributions with the IDL function 
CONVOLVE. Locations on the resulting image correspond- 
ing to negative extinction are set to zero. The point of 
peak signal provides the shift required to best align the 
isochrone model to the data. 

Although this is not a statistically rigorous method for 
isochrone fitting, this cross-correlation method is concep- 
tually straightforward, simple to code, it can efficiently 
test hundreds of models in an automated fashion to more 
quickly cover the age - composition parameter space (a 

32 http:/ /stev. oapd.inaf.it/cgi-bin/cmd 
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Figure 15. Left: T e g- - log<? diagram of five stars with Keck/HIRES spectra and properties derived with SME. Numbers indicate CWW 
ID. Padova isochrones overlaid with [M/H] = +0.064 and log* = 9.3 (blue, 2 Gyr), 9.4 (red, 2.51 Gyr), and 9.5 (green, 3.16 Gyr). 
Isochrones with ages of 2.34 and 2.7 Gyr encompass the error bars of CWW 44. The red isochrone at 2.51 Gyr shows the best fit to these 
spectroscopic properties; we also derive this age solution in our isochrone fits to NIR and optical photometric color-magnitude diagrams. 
Right: CFHT/MegaCam (g' — i') CMD for the same five stars, wit h 0.0 3 mag. error bars. Padova isochrones of same age and color scheme 
are overlaid with m — M = 7.35, Ay = 0.25. For a discussion, see £14.31 
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Figure 16. Inferring distance and extinction: brute force least squares fit of g'r'i'JHKs photometry for CWW 91 to a 1.03 Mq star 
with SED drawn from a Padova isochrone with logt = 9.4 and [M/H] = +0.064 (mass, age and metallicity suggested by SME analysis 
of a Keck/HIRES spectrum). We calculate \ 2 = S (x< — /& 2 , where £; is the measurement, fj,i is the model value, and assuming 
a g'r'i' = 0.03 mag. and ctjhk s = 0.025 mag. \ 2 i g minimized at m — M = 7.44 ± 0.02 and Ay = 0.25 ± 0.03. See i|4,3l for a discussion. 



few minutes), and it provides a diagnostic for model se- 
lection: the model with the maximum cross-correlation 
signal. In the next section, we will demonstrate that 
this method provides solutions essentially identical to the 
r 2 maximum likelihood method. Three panels in Figure 
[171 plot age, distance and extinction versus the cross- 
correlation signal, normalized to the maximum value, 
with metallicity fixed at [M/H] = +0.064. The results 
from fitting Dartmouth models arc also plotted, and 
offset by +0.06 for clarity: see Although broad, 

there are clear peaks in each diagram at \ogt = 9.4 
(2.51 Gyr), m — M = 7.32 (291 pc), and A v = 0.23 



or E(B — V) = 0.075 assuming R v = 3.1. 

The fourth panel of Figure [T7J plots the isochrones 
(gray) of each solution for logt = 9.1 to 9.56, in steps 
of 0.01. The R147 members used in the fit are plotted in 
black. The best model quoted above is overlaid in yellow, 
and two additional models are also included at younger 
and older ages, providing points of reference for how the 
fits appear to rotate across the CMD counter clockwise 
with increasing age. This rotation is primarily caused 
by the overabundance of cluster members at the MSTO, 
which serves to anchor each fit to the MSTO location. 
Those models which also pass through the handful of K 
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dwarfs and red giant branch are rewarded with a higher 
cross-correlation signal, creating the broad peaks in the 
other panels of Figure [17] Perhaps if the K dwarfs and 
red giants were weighted more heavily than the turnoff 
stars we would see more strongly peaked results. 

We ran the optical fit with [M/H] = +0.08 and find 
the best model is log* = 9.39, m — M = 7.35, Ay = 
0.26. The optical results are presented in Figure [TBI The 
consistency between the optical and NIR CMD fits and 
the parameters obtained with the T e g - logg diagram 
further justifies our use of these prelimnary optical data. 

4.5. Validating our method with r 2 

iNavlor fc Jeffriesl ()2006l ) have developed a maximum- 
likelihood method called r 2 for fitting model isochro nes 
to color - magnitude diagrams (see also Naylor 2009) 33 . 
This method simulates a cluster CMD from an isochrone, 
with a user-defined binary fraction. The user supplies a 
star catalog including the color, magnitude, and photo- 
metric errors which the t 2 code assumes are normally 
distributed. This method is powerful because it natu- 
rally accomodates errors in both color and magnitude, 
and accounts for the binary sequence. The code per- 
forms a grid search across a specified range of distances 
and ages, for isochrones of a given metallicity and red- 
dening, and identifies best values, confidence intervals, 
and returns two diagnostics for assessing how well the 
model describes the data: a reduced-r 2 (analogous to 
xt) and a probability value, Pr. 

While the r 2 code includes an isochrone library, we 
make use of the user-supplied isochrone feature and pass 
it the Padova grids, with [M/H] = +0.064. The r 2 code 
does not currently solve for Ay , so we de-redden our cat- 
alog before running r 2 , and iterate for a range of redden- 
ing values. We selected the 56 R147 members of highest 
confidence ('Y'), excluding the later K dwarfs and red 
giants. We ran r 2 for Ay = O.Oto 0.5, and allowed r 2 
to search distances ranging from 200 to 400 pc (range 
suggested by the HIP1 (-270 pc) and HIP2 (~200 pc) 
parallaxes, plus a little extra on the far side), and ages 1 
to 4 Gyr (step size is 0.01 in logi, encompassing the 2.5 
Gyr value suggested by our fit to the T e g - log g diagram 
in Figure IT~5]) . 

We find high probabilities for a large suite of models 
demonstrating the flatness of the r 2 space and high de- 
gree of degeneracy between age, metallicity, distance and 
extinction. The degeneracy is accentuated in this partic- 
ular case because we had to remove the saturated red 
giant branch from the fit. Figure [19] plots the best dis- 
tance and age values for the range of extinctions. If we 
apply the age constraint from the T e g - log g diagram fit 
(2.25 to 2.7 Gyr, illustrated by the black error bar on the 
right side of Figure fl9]> . then this restricts distance and 
extinction to m - M = 7.25 - 7.4 and Ay = 0.2 - 0.3. At 
2.5 Gyr, Ay = 0.25 and m—M = 7.35, the r 2 maximum- 
likelihood method derives values equal to those we deter- 
mined with our 2D cross-correlation method. 

4.6. Comparision between Padova, PARSEC and 
Dartmouth isochrone models 



Many groups are developing stellar evolution models. 
While there is remarkable agreement between these mod- 
els, there are noticeable departures especially around 
the turnoff and red giant branch due to choice of in - 
put physics and Solar composition ([Potter et al.l [2008). 
In fact, the model choice is the greatest source of uncer- 
tainty when determining the fundamental properties of 
star clusters via isochrone fitting. 

When we began our analysis, the Padova group pro- 
vided the only models with synthetic MegaCam pho- 
tometry. Since then, Dartmouth models have become 
available. Most recently, Padova has issued updated 
isochrones, which are now referred to as PARSEC models 
(|Bressan et al.|[20T2h : they have revised the major input 
physics, lowered the Solar metallicity from Z = 0.019 
to Z ~ 0.015, and now include the pre-main sequence 
(irrelevant for this work). 

The Padova (now PARSEC) webtool parameterizes 
composition with Z, where [M/H] = log(Z/Z©). We 
used Z = 0.022 and 0.017 to query Padova and PARSEC 
models at [M/H] w +0.065. The Dartmouth webtool 34 
uses [Fe/H] instead of Z, and we set [Fe/H] = +0.1, ac- 
cording to the SME result for CWW 91. The left panel 
of Figure l20l shows that isochrones from Padova at 2.51 
Gyr, PARSEC at 3.25 Gyr and Dartmouth at 3 Gyr map 
out the same sequence on the T e e - log g diagram, except 
for the low mass departure at ~0.8 M Q already noted. 
The middle and right panels plot the optical and NIR 
isochrones, with Ay = 0.25 and m — M = 7.35. The dif- 
ferences between models introduces an additional age un- 
certainty of at least 750 Myr. Assuming PARSEC more 
accurately models stellar evolution than their Padova 
predecessor, the difference between PARSEC and Dart- 
mouth reduces the age uncertainty to only — 250 Myr. 
Figure [17] displays the results from fitting the 2MASS 
(J - K s ) CMD with Dartmouth [Fe/H = +0.1 models, 
illustrated with the square symbols, and offset vertically 
from the Padova results by +0.06 for clarity. The peak 
is shifted to older ages relative to Padova, while the dis- 
tance and extinction remain basically consistent - a con- 
sequence of the similarity between the 2.5 Gyr Padova 
model and 3 Gyr Dartmouth model in T e g - logg space, 
and each group's color - temperature transformations. 

We will explore additional models (starting with 
BaSTI: Bag of Stellar Tracks and Isochrones, available 
at http://albione.oa-teramo.inaf.it/) in a more detailed 
analysis in a future work, where we intend to perform 
a simultaneous 7-band isochrone fit using our optical 
g'r'i'z' 2MASS JHK S , and our UKIRT JK photom- 
etry. 

4.7. Distance: photometric versus astrometric 

IDias et al.1 pOOl identified HIP 94635 (CWW 1) and 
HIP 94803 (CWW 2) as kinematic members of Ruprecht 
147, and used their HIP1 parallaxes to derive a distance 
of 270 pc to the cluster. HIP1 lists tt = 3.57 ± 1.01 
and 3.75 ± 1.04 mas for these stars respectively, while 
HIP2 provides larger values of 5.48 + 0.65 and 4.9 2 + 0.79. 
Disregarding the advice of Ivan Leeuwenl (|2007bT ) against 
deriving distances and distance moduli from parallaxes 
when the relative error is greater than 10%, one would 



code available at http:/ /www. astro. ex. ac.uk/people/timn/tau-squared/'' http:/ /stellar. dartmouth.edu/models/ 
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m - M (mag) J - K 

Figure 17. Isochrone fitting to NIR CMD with our 2D cross-correlation method: We tested Padova models (circles) with logi = 9.1 to 
9.6 (0.01 step size, Ril.25 to 4 Gyr) and [M/H] = +0.064, and Dartmouth models (squares) with [Fe/H] = +0.1 with ages running from 2 
to 3.5 Gyr. Our technique simulates stars clusters with 10 4 stars, including binaries, and computes the m — M and Ay that best matches 
the R147 photometry. The model with maximum signal in each age / visual extinction / distance modulus bin (0.01, 0.05 mag., 0.025 
mag.) is plotted in each panel, color coded according to age. We find that the model with the maximum cross-correlation signal has an age 
of 2.45 Gyr, m — M = 7.34 and Ay = 0.26. Bottom — Right: All solutions arc plotted in gray along with R147 photom etry as black 
points. Three models are highlighted: the best fit model (yellow), along with an older and younger model for reference. See i|4.4l for more 
discussion. 



determine a distance from HIP1 of 270 pc, and 193 pc 
from HIP2. 

We also identify HIP 94435 (CWW 13) as a member 
of R147. Although HIP1 lists a discrepant parallax ir = 
2.42 ± 1.25 (413 pc), HIP2 gives 4.64 ± 1.19 (216 pc) 
consistent with the other 2 HIP stars. While the HIP2 
results are self-consistent, Figure [2T1 demonstrates there 
is simply no way an isochrone can be drawn through 
either the optical (g' - i') or NIR (J - K s ) CMDs at the 
HIP2 distance moduli. The HIP2 parallax distances are 
all too close by wlOO pc, compared to the photometric 
distances. This is reminiscent of the Pleiades problem, 
where the HIP1 parallax placed the cluster 0.23 mag. 
closer than distances inferred from main sequence fitting 



and other methods. iSoderblom et al.1 (|2005| ) utilized the 
HST Fine Guidance Sensor to measure a new parallax 
distance consistent with the other non-HIP results. 

5. FINAL SYNTHESIS AND SOURCES OF UNCERTAINTY 

We have fit Padova isochrone models to three sep- 
arate datasets: a T c r - logg diagram consisting of 5 
stars with values derived from SME, and both NIR 
(J - Ks) and optical {g 1 - i') CMDs. We have cho- 
sen to use the Padova isochrones for our preliminary 
investigation into the propeties of R147 because they 
provide colors in the MegaCam filter set. Other mod- 
els (e.g. Dartmouth and Yonsei-Yale) show differences 
in the main sequence turnoff region, which provides the 
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Figure 18. Similar to Figure [T71 except we fit the optical CFHT/MegaCam (</' — i') color - magnitude diagram, with the saturated red 
giants discarded (open circles). The lower K dwarfs were also rem oved (also open circles), because the Padova models are known to run 
hotter than the main sequence at ~0.8 Mq, as illustrated in Figure 1201 We find that the model with the maximum cross-correlation signal 
has an age of 2.45 Gyr, m — M = 7.35 and Ay = 0.26. The isochrones plotted in the bottom - right panel have the following properties: 
Red - log* = 9.4, [M/H] = +0.08, m-M = 7.35 and Ay = 0.25. Blue - logt = 9.53, [M/H] = +0.08, m-M = 6.92 and Ay = 0.11. See 
£14.41 for more discussion. The fact that we derive essentially identical properties from our spectroscopi c T r ff - logg diagram, and our NIR 
and optical CMDs validates the accuracy of the CFHT Elixir and TERAFTX photometric reduction ( 8)2.5^ . 



primary age constraint; and in the lower main sequence, 
where Padova runs bluer than Dartmouth and Yonsei- 
Yale. Our results, especially for age, therefore depend 
heavily on our chosen model. 

Our solution is subject to additional sources of uncer- 
tainty, including photometric error, unresolved multiple 
star systems, and the possibility of differential extinc- 
tion. Main sequence fitting, both to the single star and 
equal mass binary sequences, provides the primary con- 
straint on the sum of distance and extinction. Ideally, 
these sequences would be vertically offset by 0.75 mag- 
nitudes (i.e. double the brightness), but for R147, there 
are separations of 1 magnitude or more. This unexpected 
offset could be explained by differential extinction, which 
would widen the sequence in both directions, by a pop- 



ulation of triple systems, or if the stars are not actually 
cluster members. In this preliminary investigation, we 
have not yet untangled this 8 Ay - multiplicity - mem- 
bership degeneracy on a star by star basis, but such an 
analysis would improve the precision of the cluster pa- 
rameters. 

At ^300 pc and 1°.25 in angular radius, we expect a 
physical radius of ~5 pc, which introduces a differential 
distance modulus of S(m — M) = 0.02 across the cluster. 
This is comprable to photometric error, and should be 
part of a more comprehensive error analysis. 

While fitting the fitting the T e s - logg diagram, as- 
suming [M/H] = +0.064, we find logt = 9.4 ± 0.03, 
oi t = 2.5 ± 0.02 Gyr. If we increase or decrease the 
metallicity by ±0.02 dcx, the age error bars remain sim- 



23 



2.6 



2.2 



\ = 0.20 



\ = 0.22 

I \ = 0.25 

H- 



\ - 0.27 

j A,, = 0.29 



A,- 0.3 



7.25 7.30 7.35 7.40 7.45 7.50 
Distance Modulus (mag.) 

Figure 19. This plot illustrates the degeneracy between extinc- 
tion, age and distance in isochrone fits to broadband photome- 
try. Results are plotted from fitting 7 different Ay values with 
t 2 , with [M/H] fixed at +0.064 according to our SME analy- 
sis. All fits returned high r 2 probabilities, and so model selec- 
tion is only possible when we place a ddit ional constraints from 
our spectroscopic analysis. In Figure 1151 we showed that a 5 
star T c g - log g diagram was best fit by a Padova isochrone with 
age = 2.5 ± 0.25 Gyr, illustrated by the error bar on the right 
side of this figure. This corresponds to a distance modulus of 
m — M = 7.35 + 0.05 — 0.1, shown by the error bar at the bot- 
tom of the figure, and Ay = 0.25 + 0.08 — 0.05. These properties 
are identical to the values we derived when fitting both the optical 
and NIR CMDs with our 2D cross-correlation techinque, validat- 
ing both our fitting method and the optical photometric reduction. 
See i|4.5l for more discussion. 



ilar and the best value for \ogt shifts by 0.02. When 
we performed the brute force SED fit to the early G 
dwarf for distance and extinction, we perturbed [M/H] 
by ±0.02 dex, log t by ±0.5, and T eff by ±50K, and found 
m-M = 7.44 ± 0.06 and A v = 0.25 ± 0.04. 

For a given metallicity and extinciton, t 2 returns typi- 
cal uncertainties of 40 - 100 Myr in age and 0.03 - 0.05 in 
distance modulus. Although the r 2 code does calculate 
two diagnostics useful for model selection, the reduced t 2 
and a probability value, for the range of parameters we 
searched, the high degree of degeneracy between the four 
cluster parameters enabled r 2 to find solutions which de- 
livered high probabilities (>60%) and reduced r 2 values 
all ~1. Instead, we will select our preferred parameter set 
by fixing the metallicity according to our spectroscopic 
SME results: [M/H] = +0.065, and the age according to 
the T e ff - \ogg result at 2.5 Gyr. This breaks the degen- 
eracy and we can then accept the corresponding distance 
and and visual extinction from r 2 as best values: m — M 
= 7.35 (d = 295 pc) and A v = 0.25. 

This is essentially identical to our 2D cross-correlation 
results, where we found peak values at an age of 2.45 
Gyr, m - M = 7.34 and A v = 0.26, for the [M/H] = 
+0.065 case. 

We set our preferred values with generous error bars at 
age = 2.5 ± 0.25 Gyr (logf = 9.4±0.03), [M/H] = 0.07 ± 
0.03, m-M = 7.25 to 7.45 (d = 280 to 310 pc), and Ay = 
0.20 to 0.30. We set the metallicity error bars according 
to our SME analysis, ignoring the one anomalously low 



metallicity result, and the hottest star with a poor xt fit. 
We set the age and error bars according to our fit to the 
T e ff - log<7 diagram (Figure [TS"]) . which is corroborated 
by the 2D cross-correlation fit (Figure fT8|). We set the 
Ay range and lower bound m — M value by placing these 
metallicity and age constraints on our r 2 results (Figure 
[ll?]) . We extend the upper bound on m-M past 7.4 
to 7.45 mag. to encompass the result of our SED fit 
to the G0/2 dwarf. The se results are summarize d in 
Table H with the value s of IKharchenko et al.l 1)20051 ) and 
iPakhomov et al.l ()2009l ) provided for comparison. 

6. SUMMARY, DISCUSSION AND UPCOMING WORK 

Over 170 years passed since Herschel first cata- 
loged Ruprecht 147 before astronomer s finally in vesti- 
gated its properti es and membership. iDias et all and 
IKharc henko et ah demonstrated that a group of 20 - 40 
stars at the location of R147 were in fact moving together 
in the plane of the sky, and estimated this group's prop- 
erties, although their analysis was hindered by (B — V) 
photometry with a li miting magnitude near the main se- 
quence turnoff. While IKharchenko et all p005h was able 
to determine an age (2.45 Gyr) from the MSTO consis- 
tent with the results of our analysis, the (B — V) main 
sequence is dominated by photometric error and there- 
fore provides a weak constraint on the distance, which 
their isochrone fitting has apparently placed 125 pc too 
close, at 175 pc compared to the 300 pc distanc e we 
find here. Neverth eless, these works by IDias et all and 
IKharc henko et ah are sig nificant because they e ssen- 
tially re-discovered R147. IPakhomov et al.l (|2009ft first 
spectroscopically determined the composition for 3 red 
giant members, showing the cluster to be super-Solar 

We queried the NOMAD catalog for stars within 1°.5 
of the cluster center, and out of the 750,000 stars, we find 
1348 with proper motion withi n 5 mas/yr of the R14 7 
value (astrometric values from IKharchenko et aLl feOOS) . 
We conducted an initial radial velocity survey at Lick 
and Palomar Observatories and for the first time con- 
firm that over 100 stars are likely members of Ruprecht 
147 and they are indeed moving together in three di- 
mensions through the Galaxy ( ^2.2.31 13~2"j) . We followed 
up this initial survey with high resolution and signal-to- 
noise Ca II H & K spectra with MMT/Hectochelle, and 
used these second epoch RVs, at higher precision, to in- 
vestigate binarity ( $2A\ 13.4.31 I3.4.4|) . We have imaged 
the cluster in four optical bands f ^2.5|) . and combined 
with 2MASS NIR photometry ( ^3.4.21) . used the result- 
ing CMDs to establish a membership list with 81 high- 
confidence members, 21 possible members, and 6 unlikely 
members ([}3|). 

We have obtained high-resolution, high-SNR spectra 
of five members f §2.3|) . and determine the metallicity to 
be super-Solar using the SME spectral synthesis code, 
and find [M/H] = +0.07 ± 0.03 and zero a-enhancement 
(» 

We have fit Padova isochrones to the T e g - log g dia- 
gram resulting from our spectroscopic analysis, and find 
that the age of R147 is best fit by a Padova isochrone 
with age of 2.5 Gyr ( N4.3)) . T c g - log g diagram isochrone 
fitting is independent of distance and visual extinction, 
which makes it a powerful tool, but this also means it 
does not directly provide any information on these pa- 
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Figure 20. A Padova isochrone (blue, 2.5 Gyr, [M/H] = +0.064) is compared to the newly released PARSEC model(green, 3.25 Gyr) and 
a Dartmouth model (red, 3 Gyr), with metallicity set at [M/H] +0.06. The Dartmouth age parameter must be increased by 500 Myr 
to match the Padova main sequence turnoff. The Padova lower main sequences are known to run hotter than many other stellar evolution 
models, diverging at approximately 0.8 Mq. While this causes the Padova optical main sequence to undershoot the K dwarfs, the NIR 
sequence well matches the Dartmouth model until 0.7 Mq. No stars below this mass were surveyed in this initial study, so the Padova 
models can be used to fit the entire NIR main sequence. 




Figure 21. The R147 optical and NIR CMDs are plotted, along with Padova isochrones: 2.5 Gyr, [M/H] = +0.064, A v = 0.25, and 
distance moduli corresponding to the HIP2 parallaxes for the 3 Hipparcos members: HIP 94635, 94803, and 94435 (CWW 1, 2, and 13). 
While the HIP2 parallaxes are all approximately in agreement, there is simply no way to place an isochrone at the distances implied by 
these parallaxes. We determine a photometric distance from opti cal and NlR isochrone fitti ng th at is about 100 pc further than the HIP2 
parallax distances, reminiscent of the Pleiades distance problem (Soderblom et al. 2005, see q4.7l l. 
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Table 4 

Table of R147 cluster properties. 



fi a [i$ RV a Age b Distance Ay c Metallicity Reference 

mas/yr mas/yr km s _1 Gyr pc mag. 

-1.1 -27.3 41.1 2.51 ± 0.25 295 ± 15 0.25 ± 0.05 [M/H] = +0.07 ± 0.03 This work 

-0.6 -27.7 41 2.45 175 0.47 ... Kharchenko et al. (20051 

40.5 1.26 ± 1.16 280 ± 100 0.34 Fei = 0.16, Fen = 0.08 Pakhomov et al. f2009V 



Note. — Our proper motions arc median values for 'Y' and 'P' members. 
a Our RV is the average of our 4 Keck velocities. Kharchenko et al.'s velocity is from Wilson (1953). Pakhomov ct al.'s velocity is 
the average of the two apparent single stars. Sec Scction |3.2| 

All three groups determined ages from Padova isochrones. Our age is the best fit parameter from a fit to a T c ^ - log g diagram 
( fl4.3t and the NIR and optical CMDs and f 34.41 . This parameter is heavily model-dependent: using PARSEC and Dartmouth 
models, we instead find 3.2 and 3 Gyr respectively (^^6). IKharchenko et all fit an optical (B — V) CMD with a Solar composition 
Padova model (Figure^ fll.2|} . Pakhomov ct al. determined composition, T e fj, and log g for 3 red giants from their high resolution 
spectroscopy, interpolated Padova models to the metallicitics for each star, and fit for mass and age using their derived T Q ff and log g 
values. 

c Ay assumes a Rv =3.1 reddening law. 
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rameters. We queried the best fit Padova isochrone for 
a star with T e g and log g closest to the values for the 
early G dwarf (CWW 91) we derived with SME, then 
performed a brute force SED fit to the resulting syn- 
thetic g'r'i' JHKg photometry, and find a minimum x 2 
at to - M = 7.44 (308 pc) and Ay = 0.25. 

We consulted the dust map of iSchleeel eFaLl (fl99l to 
set an upper limit on the amount of visual extinction to- 
ward R147, A v < 0.5 (SJHD, tnen nt Padova isochrones 
to the (g' — i') and (J — Ks) CMDs using a 2D cross- 
correlation technique devel oped he r e, whi ch was inspired 
by and validated with the INavlorl (|2009D t 2 maximum- 
likelihood method f ^4.4j) . We find that without addi- 
tional constraints from spectroscopy or additional pho- 
tometric bands, just fitting a single CMD with isochrone 
models yields a suite of solutions all with high r 2 prob- 
abilities, due to the high degree of degeneracy between 
age, composition, distance and visual extinction. If we 
break this degeneracy with the spectroscopic metallicity 
and age, we find to — M = 7.33 and Ay = 0.23 from r 2 . 

We found that the best models derived from fitting 
both (J — Ks) and (g' — i') CMDs with this technique 
corroborates the age we determined from the T e e - log g 
fit, as well as the distance and extinction corresponding 
to this age in our r 2 fits. 

We recognize significant uncertainty in our solution 
from the unresolved binary population and possibility 
of differential extinction across this large cluster. These 
results are also heavily mo del- dependent ( £ )4.6|) . 

6.1. Discussion and Future Work 

The R147 single star main sequence is not well defined, 
but blends smoothly into what is apparently the binary 
population. Evidently, R147 has a large binary fraction. 
The stellar population has encountered approximately 3 
Gyr of Galactic gravitational tidal forces. Evaporation of 
the lightest- mass members should proceed first, both be- 
cause low mass stars are preferentially ejected in 3-body 
encounters, and mass segregation gives them a larger ef- 
fective radius and more susceptibility to Galactic tides. 
This process also preferentially ejects single stars from 
the cluster, as multiple star systems of similar spectral 
type have a greater bound mass. The ongoing dynamical 
evolution of an open cluster thus tends to increase the 
binary fraction, and we may be seeing this effect in the 
R147 main sequence. 

Our membership list is top heavy, dominated by F 
stars (dwarfs, MSTO and subgiants) and red giants, with 
fewer numbers of G dwarfs, and only a handful of early K 
dwarfs. This dynamical evolution and evaporation could 
also explain the paucity of low mass members, and per- 
haps if any are left, they exist predominantly in mul- 
tiples. But this can also be explained by observational 
bias: the NOMAD and UCAC-3 proper motion errors 
increase for the fainter members. These stars also be- 
gin to blend into the Galactic background in the color 
- magnitude diagrams, further complicating candidate 
identification. We are currently conducting a radial ve- 
locity survey for lower mass members, but this question 
will only finally be settled by deriving precise proper mo- 
tions for the faint stars in the R147 field, which we in- 
tend to do by re-imaging the cluster with MegaCam in 
the near future. 

If we are able to identify single M dwarfs, these will be 



the only old (>1 Gyr), single cool dwarfs with known ages 
and compositions bright enough to admit close spectro- 
scopic study. Once the white dwarf population is iden- 
tified, it will provide an independent age estimate for 
the cluster and inform studies of white dwarf cooling 
curves. At 300 pc, chromospheric activity diagnostics 
are measureable, as is Lx, and R147 should prove use- 
ful for studying the evolution of angular momentum and 
magnetic activity at intermediate ages. 

For these reasons and more, we will continue our efforts 
to characterize Ruprecht 147 and establish it as a new 
and important benchmark for stellar astrophysics. 



J.L.C. acknowledges support from the National Science 
Foundation Graduate Research Fellowship Program sup- 
ported by NSF grant No. DGE1255832, the Stephen B. 
Brumbach graduate fellowship, and the Zaccheus Daniel 
travel grant program. AW. acknowledges a Hunter R. 
Rawlings III Cornell Presidential Research Scholarship. 

J.T.W. conceived of and oversaw the project, and col- 
lected all of the candidate stellar spectra. A.W. and 
J.T.W. conducted the radial velocity survey at Palomar 
and Lick and produced the preliminary R147 member- 
ship list. J.M.B. conducted the SME analysis. J.A.J, was 
PI on the program that obtained the optical photome- 
try from CFHT/MegaCam, and also conducted prelim- 
inary SME analysis of Lick and Keck spectra to deter- 
mine the composition of R147 (results are not reported in 
this work, though they informed this study). J.L.C. and 
J.T.W. performed the Hectochelle RV survey. J.L.C. was 
primarily responsible for the writing of the manuscript 
and the final analysis and synthesis of all data. 

We would like to thank Gabor Fiiresz and Andrew 
Szentgyorgyi for assisting the Hectochelle data reduc- 
tion; Tim Naylor for assisting us with his r 2 code; Aaron 
Dotter for providing access to the Dartmouth isochrones 
with the CFHT/MegaCam filters; David Monet and 
Stephen Levine for providing the NOMAD catalog on 
HDD; James Graham and James Lloyd for supporting 
this research; Matthew Muterspaugh for sharing and 
swapping Lick 3-m time; and Kevin Covey for providing 
extensive comments on an early draft of this paper. We 
would also like to thank all observing staff and telescope 
operators at MMT/FLWO, CFHT, Lick, Palomar, and 
Keck; and Geoff Marcy, Andrew Howard and the Cali- 
fornia Planet Survey observing team for their assistance 
in acquiring the Keck spectra. We also appreciate our 
anonymous referee for offering thorough and construc- 
tive suggestions that improved this paper. Finally, we 
thank Debra Fischer, Jeff Valenti, Adam Kraus, Steve 
Saar, S0ren Meibom, Andrew West, Kevin Covey, Mar- 
cel Agiieros, Suzanne Hawley, Ivan King, Jay Anderson, 
Bob Mathieu, Ken Janes, and Eric Mamajek for helpful 
conversations, suggestions, and support. 

This work is based on observations obtained with 
MegaCam, a joint project of CFHT and CEA/DAPNIA, 
at the Canada-France-Hawaii Telescope (CFHT) which 
is operated by the National Research Council (NRC) of 
Canada, the Institute National des Sciences de l'Univers 
of the Centre National de la Recherche Scientifique of 
France, and the University of Hawaii. Observing time 
was granted by the University of Hawaii Institute for 
Astronomy TAC. These data were reduced at the TER- 



27 



APIX data center located at the Institut d'Astrophysique 
de Paris. 

Observations reported here were obtained at the MMT 
Observatory, a joint facility of the Smithsonian Insti- 
tution and the University of Arizona. MMT telescope 
time was granted by NOAO (Project PA-10A-0378, PI 
J. Wright), through the Telescope System Instrumenta- 
tion Program (TSIP). TSIP is funded by NSF. 

This publication makes use of data products from the 
Two Micron All Sky Survey, which is a joint project of 
the University of Massachusetts and the Infrared Pro- 
cessing and Analysis Center/California Institute of Tech- 
nology, funded by NASA and the NSF. 

Some of the data presented herein were obtained at the 
W.M. Keck Observatory, which is operated as a scientific 
partnership among the California Institute of Technol- 
ogy, the University of California and the National Aero- 
nautics and Space Administration. The Observatory was 
made possible by the generous financial support of the 
W.M. Keck Foundation. The authors wish to recognize 
and acknowledge the very significant cultural role and 
reverence that the summit of Mauna Kea has always 
had within the indigenous Hawaiian community. We are 
most fortunate to have the opportunity to conduct ob- 
servations from this mountain. 

This research made use of Montage, funded by NASA's 
Earth Science Technology Office, Computation Tech- 
nologies Project, under Cooperative Agreement Number 
NCC5-626 between NASA and the California Institute of 
Technology. Montage is maintained by the NASA/IPAC 
Infrared Science Archive. 

This research made use of the WEBDA database oper- 
ated at the Institute for Astronomy of the University of 
Vienna, NASAs Astrophysics Data System Bibliographic 
Services, and the SIMBAD database and the VizieR cat- 
alogue access tool operated at CDS, Strasbourg, France. 

Any opinions, findings, and conclusions or recommen- 
dations expressed in this material are those of the authors 
and do not necessarily reflect the views of the National 
Science Foundation or any other institute cited above. 

Facilities: Shane, Hale, MMT, CFHT, Keck:I 

REFERENCES 



Agiicros, M. A. ct al. 2011, ApJ, 740, 110 

Alter, G., Ruprecht, J., & Vanysek, J. 1958, Catalogue of star 
clusters and associations, ed. G. Alter, J. Ruprecht, & 
J. Vanysek 

Archinal, B. A., & Hynes, S. J. 2003, Star clusters 
Bertin, E. 2006, in Astronomical Society of the Pacific Conference 
Series, Vol. 351, Astronomical Data Analysis Software and 
Systems XV, ed. C. Gabriel, C. Arviset, D. Ponz, & S. Enrique, 
112 

Bertin, E., Mcllicr, Y., Radovich, M., Missonnicr, G., Didelon, P., 
& Morin, B. 2002, in Astronomical Society of the Pacific 
Conference Series, Vol. 281, Astronomical Data Analysis 
Software and Systems XI, ed. D. A. Bohlcnder, D. Durand, & 
T. H. Handley, 228 

Brcssan, A., Marigo, P., Girardi, L., Salasnich, B., Dal Cero, C, 
Rubele, S., & Nanni, A. 2012, ArXiv e-prints 

Burnham, R. 1966, Burnham's celestial handbook; an observer's 
guide to the universe beyond the solar system. A descriptive 
catalog and reference handbook of deep-sky wonders for the 
observer, student, research worker, amateur or professional 
astronomer., ed. Burnham, R. 

Cardelli, J. A., Clayton, G. C, & Mathis, J. S. 1989, ApJ, 345, 
245 



Chubak, C, Marcy, G., Fischer, D. A., Howard, A. W., Isaacson, 
H., Johnson, J. A., & Wright, J. T. 2012, ArXiv c-prints 

Dias, W. S., Alessi, B. S., Moitinho, A., & Lepine, J. R. D. 2002, 
A&A, 389, 871 

Dias, W. S., Assafin, M., Florio, V., Alessi, B. S., & Li'bero, V. 

2006, A&A, 446, 949 
Dias, W. S., Lepine, J. R. D., & Alessi, B. S. 2001, A&A, 376, 441 
Dotter, A., Chaboyer, B., Jevremovic, D., Rostov, V., Baron, E., 

& Ferguson, J. W. 2008, ApJS, 178, 89 
Dreyer, J. L. E. 1888, MmRAS, 49, 1 
Drimmel, R., & Spergel, D. N. 2001, ApJ, 556, 181 
Elsanhoury, W. H., Hamdy, M. A., Nouh, M. I., Saad, A. S., & 

Saad, S. M. 2011, ISRN Astronomy and Astrophysics, 

vol. 2011, id.#127030, 2011 
Fabricant, D. et al. 2004, in Presented at the Society of 

Photo-Optical Instrumentation Engineers (SPIE) Conference, 

Vol. 5492, Society of Photo-Optical Instrumentation Engineers 

(SPIE) Conference Series, ed. A. F. M. Moorwood & M. Iye, 

767-778 

Fuhrmann, K., Pfciffcr, M., Frank, C, Rcctz, J., & Gehren, T. 

1997, A&A, 323, 909 
Furesz, G., Szentgyorgyi, A. H., & Mcibom, S. 2008, in Precision 

Spectroscopy in Astrophysics, ed. N. C. Santos, L. Pasquini, 

A. C. M. Correia, & M. Romaniello , 287-290 
Girardi, L., Bressan, A., Bertelli, G., & Chiosi, C. 2000, A&AS, 

141, 371 

Hcrschel, J. F. W. 1833, Philosophical Transactions of the Royal 

Society of London, 123, 359 
— . 1863, Royal Society of London Proceedings Series I, 13, 1 
H0g, E. et al. 2000, A&A, 355, L27 

Hora, J. L., Luppino, G. A., & Hodapp, K.-W. 1994, in Society of 
Photo-Optical Instrumentation Engineers (SPIE) Conference 
Series, Vol. 2198, Society of Photo-Optical Instrumentation 
Engineers (SPIE) Conference Series, ed. D. L. Crawford & 
E. R. Craine, 498-503 
Jones, D. O., West, A. A., & Foster, J. B. 2011, AJ, 142, 44 
Kharchenko, N. V. 2001, Kinematika i Fizika Nebesnykh Tel, 17, 
409 

Kharchenko, N. V., Piskunov, A. E., Roser, S., Schilbach, E., & 

Scholz, R. 2005, A&A, 438, 1163 
Lallement, R., Welsh, B. Y., Vergely, J. L., Crifo, F., & Sfeir, D. 

2003, A&A, 411, 447 
Landsman, W., Aparicio, J., Bergeron, P., Di Stefano, R., & 

Stecher, T. P. 1997, ApJ, 481, L93 
Lynga, G., & Palous, J. 1987, A&A, 188, 35 
Magnier, E. A., & Cuillandre, J.-C. 2002, in Society of 

Photo-Optical Instrumentation Engineers (SPIE) Conference 

Series, Vol. 4844, Society of Photo-Optical Instrumentation 

Engineers (SPIE) Conference Series, ed. P. J. Quinn, 343-350 
Marigo, P., Girardi, L., Bressan, A., Groenewegen, M. A. T., 

Silva, L., & Granato, G. L. 2008, A&A, 482, 883 
Mathieu, R. D. 1983, PhD thesis, California Univ., Berkeley. 
Mathicu, R. D., Latham, D. W., & Griffin, R. F. 1990, AJ, 100, 

1859 

Mcibom, S. et al. 2011, ApJ, 733, L9 

Mcrmilliod, J. 1995, in Astrophysics and Space Science Library, 
Vol. 203, Information & On-Line Data in Astronomy, ed. 
D. Egret & M. A. Albrecht, 127-138 

Mermilliod, J., & Paunzen, E. 2003, A&A, 410, 511 

Mink, D. J., Wyatt, W. F., Caldwell, N., Conroy, M. A., Furesz, 
G., & Tokarz, S. P. 2007, in Astronomical Society of the Pacific 
Conference Series, Vol. 376, Astronomical Data Analysis 
Software and Systems XVI, ed. R. A. Shaw, F. Hill, & 
D. J. Bell, 249-+ 

Monet, D. G. et al. 2003, AJ, 125, 984 

Naylor, T. 2009, MNRAS, 399, 432 

Naylor, T., & Jeffries, R. D. 2006, MNRAS, 373, 1251 

Nidever, D. L., Marcy, G. W., Butler, R. P., Fischer, D. A., & 
Vogt, S. S. 2002, ApJS, 141, 503 

Norris, J. M., Wright, J. T., Wade, R. A., Mahadevan, S., & 
Gcttel, S. 2011, ApJ, 743, 88 

Pakhomov, Y. V., Antipova, L. I., Boyarchuk, A. A., Bizyacv, 
D. V., Zhao, G., & Liang, Y. 2009, Astronomy Reports, 53, 660 

Peri, M. L. 1995, PhD thesis, AA(California Inst, of Tech.) 



Perryman, M. A. C, & ESA. 1997, The HIPPARCOS and 
TYCHO catalogues. Astrometric and photometric star 
catalogues derived from the ESA HIPPARCOS Space 
Astrometry Mission (The Hipparcos and Tycho 
catalogues. Astrometric and photometric star catalogues 
derived from the ESA Hipparcos Space Astrometry Mission, 
Publisher: Noordwijk, Netherlands: ESA Publications Division, 
1997, Series: ESA SP Series vol no: 1200, ISBN: 9290923997 
(set)) 

Pichardo, P>., Moreno, E., Allen, C, Bedin, L. R., Bellini, A., & 

Pasquini, L. 2012, AJ, 143, 73 
Roeser, S., Demleitncr, M., & Schilbach, E. 2010, AJ, 139, 2440 
Ruprecht, J. 1966, Bulletin of the Astronomical Institutes of 

Czechoslovakia, 17, 33 
Sanders, W. L. 1971, A&A, 14, 226 

Schilbach, E., Kharchenko, N. V., Piskunov, A. E., Roser, S., & 

Scholz, R.-D. 2006, A&A, 456, 523 
Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 

525 

Skrutskie, M. F. et al. 2006, AJ, 131, 1163 

Soderblom, D. R. 2010, ARA&A, 48, 581 

Soderblom, D. R., Nelan, E., Benedict, G. F., McArthur, B., 
Ramirez, I., Spiesman, W., & Jones, B. F. 2005, AJ, 129, 1616 

Szcntgyorgyi, A. H., Cheimets, P., Eng, R., Fabricant, D. C, 
Geary, J. C, Hartmann, L., Pieri, M. R., & Roll, J. B. 1998, in 
Society of Photo-Optical Instrumentation Engineers (SPIE) 
Conference Series, Vol. 3355, Society of Photo-Optical 
Instrumentation Engineers (SPIE) Conference Series, ed. 
S. D'Odorico, 242-252 

Thackeray, A. D. 1939, MNRAS, 99, 492 



Trumpler, R. J. 1930, Lick Observatory Bulletin, 14, 154 
Valcnti, J. A., & Fischer, D. A. 2005, ApJS, 159, 141 
Valenti, J. A., & Piskunov, N. 1996, A&AS, 118, 595 
van Leeuwen, F., ed. 2007a, Astrophysics and Space Science 

Library, Vol. 350, Hipparcos, the New Reduction of the Raw 

Data 

— . 2007b, Astrophysics and Space Science Library, Vol. 350, 

Hipparcos, the New Reduction of the Raw Data 
Vogt, S. S. 1987, PASP, 99, 1214 

Vogt, S. S. et al. 1994, in Proc. SPIE Instrumentation in 

Astronomy VIII, David L. Crawford; Eric R. Craine; Eds., 

Volume 2198, p. 362, 362-+ 
Wilson, O. C, & Vainu Bappu, M. K. 1957, ApJ, 125, 661 
Wilson, R. E. 1953, Carnegie Institute Washington 

D.C. Publication, 
Wright, J. T., Marcy, G. W., Butler, R. P., & Vogt, S. S. 2004, 

ApJS, 152, 261 

Zacharias, N. et al. 2009, VizieR Online Data Catalog, 1315, 
Zacharias, N., Monet, D. G., Levine, S. E., Urban, S. E., Gaumc, 

R., & Wycoff, G. L. 2004a, in Bulletin of the American 

Astronomical Society, Vol. 36, Bulletin of the American 

Astronomical Society, 1418 — h 
Zacharias, N., Urban, S. E., Zacharias, M. I., Wycoff, G. L., Hall, 

D. M., Monet, D. G., & Rafferty, T. J. 2004b, AJ, 127, 3043 
Zhu, Z. 2009, Research in Astronomy and Astrophysics, 9, 1285 
Zombeck, M. 2007, Handbook of Space Astronomy and 

Astrophysics: Third Edition, ed. Zombeck, M. (Cambridge 

University Press) 



Table 5 

Membership list 



CWW ID 

1 

2 

3 

4 

5 

6 

7 

8 

9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 
44 
45 
46 
47 
48 
49 
50 
51 
52 
53 
54 
55 



2MASS ID 


Mra(°>) 




ec 


M 


g' 


g'-i' 


J a 


J-K s a 


RVlp 


RV H 


Mem. b 


MemFlag c 


Notes d 


19152612-1605571 


-1.0 


(0.7) 


-27 


4 


(o 


6 ) 


8.37 


1.18 


5.31 


0.76 


38 


5 


- 


Y 


YP-YY 


RG 


19172384-1604243 


-2.2 


(0.8) 


-27 


6 


(0 


6 ) 


8.66 


1.26 


5.27 


0.83 


43 


4 


- 


Y 


YY-YY 


RG 


19161966-1634094 


-1.5 


(1.6) 


-25 


4 


(1 


5) 


8.65 


-0.10 


8.05 


-0.02 


- 




- 


Y 


Y- -BB 


BS 


19171130-1603082 


-0.9 


(1.5) 


-29 


1 


(1 


5) 


9.20 


1.10 


6.40 


0.77 


42 


7 


41.1 


Y 


YYYYY 


RG 


19164073-1616411 


-0.3 


(1.6) 


-25 


4 


(1 


5) 


8.81 


-0.04 


7.95 


0.05 


- 




- 


Y 


Y- -BB 


BS 


19170343-1703138 


-0.7 


(1.9) 


-30 


1 


(1 


8 ) 


9.21 


1.02 


6.42 


0.71 


46 


2 


- 


Y 


YY-YY 


RG 


19183747-1712575 


-0.4 


(1.9) 


-26 


7 


(1 


8 ) 


- 


- 


6.41 


0.67 


42 


4 


- 


Y 


YY-Y- 




19181439-1641226 


1.8 


(1.5) 


-25 


3 


(1 


5) 


8.93 


0.33 


7.66 


0.18 


- 




- 


Y 


Y- -BB 


BS 


19140272-1554055 


-0.4 


(1.5) 


-26 


1 


(1 


4) 


9.30 


1.09 


6.31 


0.75 


42 


1 


- 


Y 


YY-YY 


RG/SB1? 


19155129-1617591 


-2.7 


(1.6) 


-26 


8 


(1 


5) 


9.17 


0.96 


6.45 


0.66 


41 


4 


40.6 


Y 


YYYYY 


RG 


19180978-1616222 


-1.4 


(1.5) 


-27 


3 


(1 


5) 


9.39 


0.98 


6.54 


0.70 


44 


2 


- 


Y 


YY-YY 


RG 


19164388-1626239 


-1.2 


(1.6) 


-24 


6 


(1 


5) 


9.60 


0.32 


8.44 


0.14 


- 




- 


Y 


Y- -BB 


BS 


19131526-1706210 


-0.1 


(1.1) 


-27 


1 


(0 


8) 


- 


- 


7.28 


0.65 


46 


4 


- 


Y 


YY-Y- 




19134817-1650059 


-1.5 


(1.4) 


-26 


2 


(1 


3) 


9.88 


0.87 


7.63 


0.56 


43 


6 


- 


Y 


YY-YY 


RG 


19164574-1635226 


-1.6 


(1.6) 


-27 


1 


(1 


5) 


9.96 


1.02 


7.38 


0.66 


46 


1 


41.4 


Y 


YYYYY 


RG 


19164823-1611522 


1.9 


(1.6) 


-26 


1 


(1 


™) 


9.96 


0.23 


9.11 


0.17 


- 




- 


Y 


Y- -BB 


BS 


19165670-1612265 


-1.7 


(1.7) 


-26 


8 


(1 


8 ) 


10.07 


0.54 


8.59 


0.37 


44 


2 


40.6 


Y 


YYYYY 




19193373-1658514 


0.3 


(4.2) 


-26 


2 


(4 


™) 


- 


- 


8.68 


0.36 


47 


2 


- 


Y 


YP-Y- 




19161456-1624071 


-0.2 


(2.0) 


-28 


6 


(2 


°) 


10.32 


0.95 


8.00 


0.62 


43 


8 


43.9 


P 


YYPYY 


RG 


19160865-1611148 


-3.4 


(1.6) 


-29 


3 


(1 


8) 


10.25 


0.57 


8.63 


0.33 


41 


8 


41.7 


Y 


YYYYY 




19132220-1645096 


-5.4 


(1.5) 


-29 


1 


(1 


6 ) 


10.28 


0.51 


8.74 


0.37 


41 


9 


- 


Y 


YY-YY 




19172382-1612488 


-1.7 


(1.7) 


-30 


8 


(1 


8) 


10.31 


0.42 


9.26 


0.44 


51 


9 


38.2 


P 


YNNYY 


SB1 


19154269-1633050 


-1.8 


(2.2) 


-30 


3 


(2 


3) 


11.67 


0.44 


10.38 


0.32 


34 


7 


41.2 


Y 


YNYYY 


SB1? 


19172865-1633313 


1.1 


(1.7) 


-27 


6 


(1 


8) 


10.29 


0.40 


9.11 


0.21 


41 


8 


- 


Y 


YY-YB 


BS? 


19133648-1548104 


-1.4 


(1.4) 


-28 


5 


(1 


5) 


10.53 


0.80 


8.34 


0.61 


40 


7 


- 


Y 


YY-YY 


RG 


19153282-1620388 


0.2 


(1.8) 


-27 


1 


(2 


0) 


10.51 


0.44 


9.03 


0.31 


46 


1 


42.0 


Y 


YYYYY 




19171984-1607383 


-2.0 


(1.9) 


-30 


8 


(2 


°) 


10.46 


0.58 


8.96 


0.29 


48 


3 


48.4 


P 


YPNYY 




19152638-1700159 


-2.1 


(2.0) 


-29 


6 


(2 


3) 


10.53 


0.56 


9.08 


0.24 


43 


9 


- 


Y 


YY-YY 




19173931-1636348 


1.2 


(1.6) 


-25 


4 


(1 


6) 


10.47 


0.57 


9.00 


0.33 


41 


6 


- 


Y 


YY-YY 




19155841-1615258 


-2.8 


(1.9) 


-28 


9 


(2 


1) 


10.57 


0.52 


9.20 


0.39 


41 


4 


40.7 


Y 


YYYPY 




19195154-1603583 


-2.4 


(2.0) 


-27 


3 


(2 


2) 


- 


- 


9.05 


0.34 


41 


7 


- 


Y 


YY-Y- 




19151540-1619517 


-1.8 


(1.8) 


-26 


8 


(2 


1) 


10.70 


0.52 


9.30 


0.33 


45 


9 


42.2 


Y 


YYYYY 




19181155-1629141 


1.4 


(1.9) 


-27 





(1 


9) 


10.66 


0.44 


9.39 


0.29 


41 


2 


- 


Y 


YY-YY 




19165477-1702129 


3.2 


(3.2) 


-28 


7 


(3 


5) 


10.63 


0.46 


9.41 


0.31 


45 


8 


- 


Y 


YY-YY 




19163976-1626316 


0.1 


(1.9) 


-24 


8 


(2 


1) 


10.53 


0.55 


9.10 


0.32 


45 


1 


41.1 


Y 


YYYYY 




19153626-1557460 


-0.7 


(1.9) 


-27 


2 


(2 


1) 


10.49 


0.53 


9.07 


0.31 


46 


4 


41.8 


Y 


YYYYY 




19163344-1607515 


-2.5 


(2.1) 


-27 


6 


(2 


2) 


10.64 


0.49 


9.32 


0.36 


34 


6 


41.0 


Y 


YNYYY 




19142651-1606340 


-2.5 


(2.0) 


-27 


2 


(2 


1) 


10.66 


0.48 


9.33 


0.41 


45 


1 


41.2 


Y 


YYYPY 




19150275-1609405 


-5.2 


(1.9) 


-28 


3 


(2 


1) 


10.66 


0.52 


9.19 


0.31 


42 


7 


41.7 


Y 


YYYYY 




19163339-1620215 


-3.0 


(1.9) 


-24 


9 


(2 


1) 


10.62 


0.56 


9.18 


0.33 


42 


1 


41.5 


Y 


YYYYY 




19170481-1636526 


4.3 


(2.0) 


-28 


6 


(2 


2) 


10.61 


0.44 


9.33 


0.30 


45 


9 


- 


Y 


YY-YY 




19183120-1614421 


-1.8 


(2.2) 


-26 


9 


(2 


2) 


10.64 


0.45 


9.32 


0.33 


44 


6 


- 


Y 


YY-YY 




19180054-1636016 


1.3 


(1.9) 


-26 


3 


(2 


0) 


10.70 


0.46 


9.44 


0.27 


44 


9 


- 


Y 


YY-YY 




19164495-1717074 


2.7 


(2.3) 


-27 


2 


(2 


5) 


10.97 


0.48 


9.74 


0.32 


42 


2 


- 


Y 


YY-YY 




19150860-1657412 


-2.0 


(2.0) 


-30 





(2 


3) 


10.82 


0.45 


9.50 


0.28 


45 


7 


- 


Y 


YY-YY 




19163525-1705075 


3.7 


(2.3) 


-28 





(2 


6) 


10.80 


0.46 


9.52 


0.30 


45 


5 


- 


Y 


YY-YY 




19131541-1616123 


-5.6 


(1.7) 


-30 





(1 


8 ) 


10.78 


0.45 


9.51 


0.32 


44 


5 




Y 


YY-YY 




19164662-1619208 


0.0 


(2.1) 


-26 


4 


(2 


2) 


10.72 


0.54 


9.34 


0.37 


46 


3 


42.3 


Y 


YYYYY 




19142907-1549056 


0.3 


(1.9) 


-25 


3 


(2 


1) 


10.76 


0.50 


9.36 


0.36 


46 


2 




Y 


YY-YY 




19162934-1645544 


-2.4 


(2.0) 


-27 


4 


(2 


3) 


10.97 


0.37 


9.80 


0.30 


47 


6 




P 


YP-YP 




19163620-1607363 


-3.0 


(2.5) 


-29 


5 


(2 


6 ) 


11.02 


0.46 


9.12 e 


0.35 


44 


9 




Y 


YY-YY 




19162169-1609510 


-0.9 


(2.1) 


-27 


7 


(2 


3) 


10.90 


0.47 


9.59 


0.34 


39 


2 


40.1 


Y 


YYYYY 




19163231-1611346 


-2.8 


(2.2) 


-29 


3 


(2 


3) 


10.93 


0.53 


9.60 


0.33 


46 


9 


45.2 


P 


YYNYP 




19165573-1603220 


-0.8 


(2.2) 


-30 


2 


(2 


4) 


11.39 


0.42 


10.10 


0.32 


48 


9 


42.8 


Y 


YPYYY 




19160452-1605313 


-0.6 


(2.1) 


-32 


1 


(2 


2) 


10.94 


0.44 


9.61 


0.30 


45 


4 


41.8 


Y 


YYYYY 





Table 5 — Continued 



CWW ID 



2MASS ID 



MDec 




g' 


g'-i' 


J a 


J-K S a 


RVlp 


RV H 


Mem. b 


MemFlag 


-25.8 


(2.0) 


- 


■ 


9.83 


0.33 


40.7 


- 


Y 


YY-Y- 


-28.6 


(2.3) 


11.23 


0.47 


9.93 


0.27 


41.9 


41.8 


Y 


YYYYY 


-28.1 


(1.7) 


11.27 


0.46 


9.95 


0.36 


44.6 




Y 


YY-PY 


-30.6 


(2.7) 


11.18 


0.48 


9.92 


0.24 


50.1 




P 


YP-PY 


-27.1 


(2.3) 


11.31 


0.49 


10.09 


0.38 


43.5 




Y 


YY-PY 


-29.1 


(2.3) 


11.18 


0.45 


9.93 


0.32 


42.2 


- 


Y 


YY-YY 


-24.9 


(2.3) 


11.25 


0.52 


9.84 


0.32 


45.0 


42.5 


Y 


"\ 7"^ 7"\. r~\ 7"\. T 

YYYYY 


-26.5 


(1.8) 


11.53 


0.52 


10.17 


0.33 


39.7 


41.4 


Y 


YYYYY 


-26.1 


(2.5) 


11.72 


0.64 


10.20 


0.40 


~ 




P 


Y-PP 


-23.7 


(3.7) 


13.88 


1.15 


11.45 


0.68 


- 


- 


Y 


Y-PY 


-27.5 


(2.3) 


11.71 


0.62 


10.15 


0.36 




- 


P 


Y— YP 


-29.1 


(2.8) 


11.84 


0.75 


10.08 


0.47 


34.7 




N 


YN-NN 


-29.8 


(3.7) 


13.57 


1.00 


11.25 


0.67 


- 


- 


Y 


Y-PY 


-30.6 


(2.3) 


11.38 


0.51 


9.99 


0.37 


42.4 


37.2 


P 


PYNPY 


-29.6 


(2.3) 


11.47 


0.43 


10.23 


0.30 


37.9 


38.9 


P 


YPPYY 


-27.0 


(2.4) 


11.93 


0.48 


10.57 


0.32 


41.4 


41.1 


Y 


"\ 7"^ r~\. r~\ r~\. r 

YYYYY 


-30.8 


(2.4) 


12.04 


0.71 


10.36 


0.40 


47.9 


46.6 


P 


YPNYP 


-26.8 


(2.4) 


11.73 


0.47 


10.43 


0.33 


42.5 


- 


Y 


YY-YY 


-28.6 


(1.8) 


11.55 


0.49 


10.23 


0.34 


42.1 


42.2 


Y 


YYYYY 


-25.5 


(3.6) 


13.20 


0.74 


11.43 


0.52 


39.5 


42.2 


Y 


*\ r~\ 7"\. r*\ 7"\. t 

YYYYY 


-77.7 


(3.8) 


12.69 


0.68 


11.13 


0.41 


43.5 


_ 


P 


NY-YY 


-25.0 


(3.2) 


12.12 


0.59 


10.56 


0.40 


50.5 


52.6 


K 


YNN Y Y 


-30.5 


(2.1) 


11.70 


0.51 


10.35 


0.31 


40.8 




Y 


YY-YY 


-28.6 


(3.1) 


12.60 


0.66 


11.00 


0.40 


42.3 


- 


Y 


"\ r~\ r~\. r~\ r~\. r 

YYYYY 


-29.5 


(2.5) 


12.13 


0.74 


10.38 


0.45 


41.5 




P 


YY-PP 


-25.5 


(2.4) 


11.77 


0.51 


10.44 


0.36 


41.2 


- 


Y 


YY-YY 


-29.6 


(2.4) 


11.97 


0.51 


10.55 


0.36 


47.9 


42.5 


Y 


YPYYY 


-29.5 


(3.3) 


12.23 


0.69 


10.64 


0.44 


42.0 


- 


Y 


YY-YY 


-26.6 


(2.9) 


12.41 


0.76 


10.61 


0.44 


46.7 


- 


Y 


YY-YY 


-27.7 


(3.1) 


12.82 


0.64 


11.23 


0.44 


42.6 


42.9 


Y 


YYYYY 


-28.7 


(2.9) 


12.94 


0.85 


10.98 


0.49 


44.7 


41.6 


Y 


YYYYY 


-26.8 


(3.0) 


12.83 


0.84 


10.91 


0.49 


45.0 


42.4 


Y 


"\ r~\ r~\. r~\ r~\. r 

YYYYY 


-42.8 


(3.7) 


13.02 


0.72 


11.28 


0.43 


47.0 


~ 


P 


PY-YY 


-31.9 


(2.5) 


12.72 


0.58 


11.07 


0.41 


47.3 




P 


YP-YY 


-31.8 


(2.9) 


12.44 


0.55 


10.95 


0.32 


42.0 


_ 


Y 


YY-YY 


-29.5 


(3.0) 


12.75 


0.61 


11.13 


0.36 


42.8 


42.2 


Y 


"\ r~\ r~\. r~\ r~\. r 

YYYYY 


-29.0 


(3.6) 


12.63 


0.74 


10.89 


0.45 


45.1 


25.2 


P 


YYNYY 


-27.6 


(2.9) 


13.02 


0.76 


11.29 


0.43 


41.8 


41.7 


Y 


YYYYY 


-27.1 


(3.2) 


13.26 


0.77 


11.43 


0.38 


43.5 


42.8 


Y 


YYYPY 


-30.1 


(2.9) 


12.78 


0.79 


10.90 


0.48 


40.6 


39.1 


P 


YYPYY 


-27.0 


(2.1) 


~ 


■ 


10.73 


0.35 


40.8 


- 


Y 


YY-Y- 


-27.9 


(3.6) 


13.03 


0.70 


11.38 


0.43 


40.9 


41.1 


Y 


"\ r~\ r~\. r~\ r~\. r 

YYYYY 


-29.5 


(2.9) 


13.52 


0.97 


11.38 


0.63 


43.2 


40.7 


Y 


YYYPY 


-29.4 


(3.0) 


13.66 


0.89 


11.52 


0.58 


44.8 


44.6 


P 


Y YJN Y Y 


-3.7 (3.7) 


14.20 


1.22 


11.67 


0.79 


42.2 




N 


NY-PY 


-31.6 


(3.7) 


15.43 


1.43 


12.52 


0.78 


46.3 


40.5 


Y 


PYYYY 


-34.1 


(3.8) 






11.85 


0.55 


44.3 




P 


PY-Y- 


-11.1 


(4.9) 


14.27 


1.15 


11.92 


0.68 


45.0 




N 


NY-YY 


-15.8 


(4.8) 






11.78 


0.65 


52.4 




N 


NN-Y- 


-58.2 


(3.7) 


14.39 


1.16 


11.87 


0.68 


50.5 




N 


NN-YY 


-31.9 


(3.7) 


15.05 


1.53 


12.06 


0.76 


49.3 


46.5 


P 


YPNYY 


-34.8 


(3.7) 


14.70 


1.26 


12.02 


0.66 


42.5 


42.1 


Y 


PYYYY 


-32.3 


(3.7) 


15.21 


1.35 


12.47 


0.71 


48.6 


42.8 


P 


NPYYY 



c Notes d 



56 19200522-1535360 1.1 

57 19170433-1623185 -0.2 

58 19172172-1535592 -0.4 

59 19151260-1705121 2.9 

60 19114731-1632485 0.8 

61 19145840-1650089 -4.5 

62 19164922-1613222 -0.5 

63 19152981-1551047 -1.8 

64 19152465-1651222 -1.5 

65 19164440-1615338 -3.6 

66 19150050-1614245 -1.5 

67 19151498-1720177 -0.1 

68 19180536-1646438 3.4 

69 19161864-1611305 -6.8 

70 19163827-1625039 1.1 

71 19154511-1623157 -0.1 

72 19165800-1614277 -2.6 

73 19160523-1652561 -3.5 

74 19150925-1552241 -0.1 

75 19161121-1621485 4.2 

76 19134334-1649109 6.6 

77 19150012-1605517 -0.4 

78 19160879-1524279 -3.8 

79 19142816-1620023 -2.6 

80 19162501-1632018 -2.0 

81 19151897-1639244 -2.1 

82 19152406-1621519 -1.3 

83 19134126-1610201 -5.3 

84 19141294-1554291 -1.8 

85 19165940-1635271 -2.8 

86 19160589-1629481 -0.3 

87 19160785-1610360 -4.5 

88 19162477-1710375 -4.3 

89 19173402-1652177 -1.6 

90 19163672-1713101 -0.1 

91 19164725-1604093 -2.2 

92 19164417-1612222 -1.5 

93 19162203-1546159 1.5 

94 19152141-1600107 -2.6 

95 19170128-1609423 -1.5 

96 19151156-1726308 -0.5 

97 19170285-1605166 -1.7 

98 19162656-1614545 0.4 

99 19161757-1600177 -2.4 

100 19145199-1541379 0.2 

101 19153354-1625368 -9.0 

102 19124958-1550340 3.4 

103 19134512-1619340 19.9 

104 19193779-1618312 11.4 

105 19181352-1614496 5.8 

106 19163680-1623032 -3.0 

107 19163732-1600050 -2.9 

108 19172940-1611577 9.1 



1.9) 
'2.2) 
1.7) 
2.5) 
2.2) 
2.2) 
2.2) 
1.7) 
2.4) 
3.7) 
2.2) 
2.5) 
3.7) 
'2.2) 
2.2) 
'2.2) 
2.3) 
2.4) 
1.8) 
3.6) 
3.8) 
'3.2) 
2.0) 
3.0) 
2.4) 
2.2) 
2.4) 
3.2) 
3.1) 
3.1) 
2.9) 
2.9) 
3.7) 
2.5) 
2.9) 
3.1) 
3.6) 
2.9) 
3.2) 
2.9) 
2.1) 
3.6) 
2.9) 
3.0) 
3.7) 
3.7) 
3.8) 
4.9) 
4.8) 
3.7) 
3.7) 
3.7) 
3.7) 



SB2 
SB2 
SB2 

SB2 



SB2 



SB1 



Table 5 — Continued 



CWWID 2MASS ID MRa(o>) ^Decfa*) g' g'-i' J a J - K s a RV L p RV H Mem. b MemFlag c Notes d 



Note. — Column Notes: (1) CWW ID - This work's star identification scheme, sorted by V magnitude. CWW — Curtis, Wolfgand and Wright. (2) 2MASS ID, also provides RA and Dec 
positions (3, 4) RA and Dec proper motions in mas/yr from PPMXL catalog (5) CFHT/McgaCam g mag. (6) g — i' mag. (7) 2MASS J mag. (8) 2MASS J — Ks mag. (9) Lick / Palomar RV 
in km s _1 (10) Hcctochelle RV in km s" 1 (11) Membership probabilities (12) Membership probabilities for each criterion (13) Notes for individual stars. Vaues in parenthesis are measurement 
errors. 

a We use 2MASS aperture photometry instead of the default PSF photometry for 18 stars, based on our analysis that these stars, and only these stars, shift position on the (J — Ks) CMD 
and that they all move toward the cluster locus. No neighbors arc resolved in our optical imaging within 5". The stars are CWW 22, 24, 27, 28, 30, 37, 38, 43, 48, 49, 57, 59, 67, 68, 90, 91, 
94, and 100. 

k Membership Probability: Y — yes, highest confidence member, P — possible / probable member, N — not likely / non-member 

C Membership Criteria: proper motion radial distance from cluster value; Lick / Palomar RV, Hectochelle RV, 2MASS (J — Ks) CMD, CFHT/McgaCam (g' — i') CMD. Confidence intervals 
defined in Table [2] A 'B' flag indicates photometry consistent with blue stragglers. A dash '-' indicates no data. 

^ Notes: BS — blue straggler, RG — red giant, SB2 — spectroscopic double line binary, SB1? — inconsistent RVs between multiple epochs 

6 Our MegaCam imaging shows CWW 51 is an optical double, with a star 1.65" away with a similar g r' i' z SED. Assuming a cluster distance of 300 pc, this angular separation translates 
into a minimum physical separation of 495 AU. This suggests that the pair actually form a wide binary, although their angular proximity could also be explained by a chance alignment. This 
double was not resolved in the 2MASS Point Source Catalog. Adding 0.75 mag. to the J band magnitude (halving the brightness, to reflect just the one star) moves CWW 51 in the (J — Ks) 
CMD to its neighbors in the (g — i ) CMD. Despite this realization, we quote the 2MASS PSC photometry here. Sec E|3.4.2| 



